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PREFACE 


This  report  presents  criteria  for  load-transfer  and  numerical  analysis 
of  axially  loaded  piles  in  clay.  These  criteria  were  selected  from  the  latest 
published  information  on  axially  loaded  piles  in  clay  and  reflect  the  current 
state-of-the-art  in  pile  foundation. 

Dr.  Andrew  G.  Heydinger,  Department  of  Civil  Engineering,  The  Univer¬ 
sity  of  Toledo,  Toledo,  Ohio,  prepared  the  report  under  Contract  No.  DACW39- 
84-M-2309.  The  work  in  developing  these  criteria  was  done  as  part  of  the  ap¬ 
plication  support  provided  by  the  Information  Technology  Laboratory  (ITL), 
formerly  Automation  Technology  Center  (ATC),  US  Army  Engineer  Waterways 
Experiment  Station  (WES),  to  the  US  Army  Engineer  Division,  Lower  Mississippi 
Valley  (LMVD).  The  point  of  contact  for  LMVD  was  Mr.  James  A.  Young. 

At  WES  the  coordination  and  monitoring  was  accomplished  by  Mr.  Reed 
Mosher,  Engineering  Application  Group,  formerly  of  Scientific  and  Engineer¬ 
ing  Application  Division  (SEAD),  ATC.  Mr.  Mosher  provided  technical  guid¬ 
ance  and  review  on  the  project  under  the  general  supervision  of  Mr.  Paul  K. 
Senter,  Acting  Chief,  Information  Research  Division  (formerly  SEAD),  and 
Dr.  W.  Radhakrishnan,  Acting  Chief,  ITL,  formerly  Chief  ATC,  WES.  Informa¬ 
tion  Products  Division,  ITL,  WES,  Editor  and  Editorial  Assistant  Mses.  Gilda 
Shurden  and  Frances  Williams,  respectively,  prepared  the  report  for 
publication. 

COL  Allen  F.  Grum,  USA,  was  the  previous  Director  of  WES.  COL  Dwayne  G. 
Lee,  CE,  is  the  present  Commander  and  Director.  Dr.  Robert  W.  Whalin  is 


Technical  Director. 


Accesion  For 

NTtS  CRA&I 
DTIC  TAB 
U.  .announced 

Justification . 


By  . — 

Di  t  ib  itiot'  | _ _ 

Avj'tability  Codes^ 
‘T  Avail  cV  d  /  or 

Di' t  jjr  cut 

U/l  i 


.5  i 


s:* 

ESs* 

V»V. / 

l.iT*  * *1  ■ 


M 


m 


»  •  *  —  »  •  *  W9  Wt 


•* 


CONTENTS 


PREPACE . 

CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC)  UNITS  OF  MEASUREMENT 
PART  I:  INTRODUCTION . 


Background . 

Purpose . 

Scope . 

PART  II:  CAPACITY  OF  PILES  IN  COHESIVE  SOILS.. 

Introduction . 

Total  Stress  Methods . 

Effective  Stress  Methods . 

Comparison  of  Predictions . 

Conclusions  and  Recommendations . 

PART  III:  PILE  SIDE  SHEAR-DISPLACEMENT  FUNCTION 

Introduction . 

Methods  for  Computing  f-z  Curves . 

Example  of  Predictive  Methods . 

Discussion  of  Predictive  Methods . 

Comparison  of  Predictive  Methods . 

Soil-Structure  Interaction  Program . 

Conclusion  and  Recommendations . 


4 

5 
5 

7 

7 

7 

17 

21 

23 

25 

25 

25 

31 

33 

40 

45 

52 


PART  IV:  ENGINEERING  PROPERTIES  OF  SOILS 

Introduction . 

Correlations  with  Plasticity . 

Normalized  Soil  Testing . 

Recommends t ions . 

REFERENCES . 

TABLES  1-6 

APPENDIX  A:  NOTATION . 


53 


53 

59 

64 

67 


A1 


CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  (met- 


ric)  units  as  follows: 

Multiply 

By 

To  Obtain 

feet 

0.3048 

metres 

inches 

0.0254 

metres 

kips  (force)  per  square  foot 

47.88026 

kilopascals 

pounds 

4.448222 

newtons 

pounds  (force)  per  square  foot 

47.88026 

pascals 

pounds  (mass)  per  cubic  foot 

16.01846 

kilograms  per  cubic  metre 

pounds  (force)  per  square  inch 

6.894757 

kilopascals 

tons  (force) 

8.896444 

kilonewtons 

i 
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RECOMMENDATIONS:  LOAD-TRANSFER 
CRITERIA  FOR  PILES  IN  CLAY 


PART  I :  INTRODUCTION 

Background 

1.  The  US  Army  Corps  of  Engineers  (CE)  designs  and/or  oversees  the  de¬ 
sign  of  many  structures  with  pile  foundations  in  cohesive  soils.  So  that 
these  structures  will  perform  satisfactorily,  not  only  the  capacity  of  the 
foundations  must  be  adequate  but  also  the  load-deformation  behavior  of  the 
foundations  must  be  predicted  to  allow  an  analysis  of  the  structures  for 
stresses  and  deformations.  The  difficult  task  of  structure  analysis  is  accom¬ 
plished  with  the  aid  of  computer  programs  which  rely  on  input  on  the  founda¬ 
tion  behavior.  The  information  in  this  report  is  intended  to  form  the  back¬ 
ground  for  developing  a  computer  subroutine  to  generate  data  to  create  design 
charts  for  piles  in  cohesive  soils  or  a  subroutine  in  a  large  interactive 
program  for  pile  design. 

2.  The  capacity  of  a  foundation  is  measured  by  its  ability  to  withstand 
loads  without  excessive  settlements.  The  capacity  of  piles  in  clay  consists 
of  the  resistance  to  penetration  at  the  tip  of  the  piles  and  the  side-shear 
resistance.  A  solution  requires  identification  of  pertinent  soils  properties 
and  the  resulting  conditions  after  the  piles  have  been  installed. 

3.  The  load-deformation  behavior  of  piles,  the  predicted  settlements, 
and  their  respective  loads  require  determination  of  functions  for  pile-soil 
interaction  behavior.  The  functions  define  the  dependence  of  side  or  tip  re¬ 
sistance  on  vertical  pile  movements.  Computer  programs  have  been  developed  at 
the  US  Army  Engineer  Waterways  Experiment  Station  (WES)  (Radhakrishnan  and 
Parker  1975)  to  make  use  of  recommended  procedures  for  predicting  pile-soil 
interaction.  However,  problems  associated  with  predicting  behavior  are  more 
complex  than  the  initial  procedures  can  account  for,  so  this  research  was 
necessary . 

4.  For  the  purpose  of  foundation  design,  soils  can  be  classified  as  co¬ 
hesive  or  cohesionless.  Cohesive  soils  are  fine-grained  soils  with  low  perme¬ 
abilities  and  shear  strengths  with  a  dependence  on  attractive  forces  between 
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the  soil  particles.  Cohesionless  soils  are  coarse-grained  soils  with  high 
permeabilities  and  shear  strengths  and  a  dependence  on  the  stresses  acting  on 
the  soil  particles.  Problems  associated  with  piles  in  cohesionless  soils  have 
recently  been  addressed  by  Mosher  (1984).  Large  pore  pressures  and  severe 
soil  remolding  occur  during  pile  installation  in  cohesive  soils,  resulting  in 
problems  unique  from  those  of  cohesionless  soils.  Such  problems  are  con¬ 
sidered  in  this  report. 


Purpose 

5.  Investigation  of  design  procedures  for  axially  loaded  single  piles 
in  cohesive  soils  and  recommendations  for  their  use  are  the  purpose  of  this 
research  project.  More  specifically,  predictive  methods  for  computing  capac¬ 
ity  and  shear  transfer  versus  vertical  pile  movement,  f-z*  curves,  were  in¬ 
vestigated.  The  specific  objective  for  this  research  was  to  investigate  piles 
of  100  ft**  or  less  in  length,  and  precast  concrete  piles  and  H-piles  that  are 
likely  to  be  used  on  CE  projects.  Additionally,  this  research  includes  infor¬ 
mation  and  recommendations  for  estimating  soil  properties  that  are  useful  for 
the  design  of  piles. 


Scope 

6.  This  report  includes  research  on  proposed  predictive  methods  for 
axial  capacity  and  on  the  effects  of  pile  installation  on  soil  properties. 

The  study  includes  the  most  recent  attempts  to  determine  the  changes  in  shear 
strength  and  the  state  of  stress  around  piles  in  clay.  The  predictive  methods 
that  are  discussed  are  total  stress  methods  based  on  in  situ  (before  pile  in¬ 
stallation)  soil  properties  and  effective  stress  methods. 

7.  Research  on  load-deformation  behavior  is  directed  toward  determin¬ 
ing  f-z  curves,  which  can  be  used  in  a  WES  program  to  compute  pile  behavior. 
Four  proposed  methods  were  analyzed  to  determine  a  recommended  method  for  CE 
purposes.  Parametric  studies  were  conducted  to  determine  the  method  that 


*  For  convenience,  symbols  and  abbreviations  are  listed  in  the  notation 
(Appendix  A). 

**  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI  (met¬ 
ric)  units  is  presented  on  page  3. 
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could  best  be  used  to  represent  actual  soil  conditions.  Results  of  four  in- 
strumented  pile-load  tests  were  also  used  to  evaluate  the  methods. 

8.  Research  on  soil  properties  consists  of  correlations  that  can  be 
used  to  determine  in  situ  soil  properties.  The  research  includes  correlations 
for  the  undrained  shear  strength,  the  effective  friction  angle,  KQ,  and 
Young's  modulus.  The  correlations  can  be  used  for  designing  piles. 


PART  II:  CAPACITY  OF  PILES  IN  COHESIVE  SOILS 


Introduction 

9-  Methods  to  compute  pile  capacity  are  referred  to  as  total  or  effec¬ 
tive  stress  methods.  Total  stress  methods  rely  on  in  situ  soil  properties  or 
stresses  and  empirical  observations  obtained  from  pile  load  tests.  Three 
methods  that  have  been  used  are  the  a-method  (American  Petroleum  Institute 
(API)  1978,  Vesic  1972),  the  8-method  (Meyerhof  1976),  and  the  X-method 
(Vijayvergiya  and  Focht  1972).  Effective  stress  methods  are  based  on  predic¬ 
tions  of  the  effective  stresses  in  the  soil  around  piles  in  clay  at  the  time 
of  loading.  The  methods  are  those  reported  by  Esrig  and  Kirby  (1979),  Kraft 
(1982),  and  Heydinger  and  O'Neill  (in  preparation).  This  section  reviews  in¬ 
formation  on  the  predictive  methods  and  determines  a  recommended  design  proce¬ 
dure  for  CE  projects. 

10.  Error  can  occur  with  these  predictive  methods  from  a  number  of 
sources.  For  piles,  the  pile  material  and  the  method  of  installation  affect 
the  results.  For  pile  testing,  there  are  differences  in  the  way  and  the  rate 
that  the  loads  are  applied.  Error  also  occurs  in  measuring  the  applied  loads, 
particularly  if  load  cells  are  not  used.  For  soils,  different  methods  have 
been  used  to  obtain  the  undrained  shear  strengths.  There  is  considerable  var¬ 
iation  in  soil  properties  for  a  given  site,  requiring  estimation  of  the  un¬ 
drained  shear  strength.  Test  results  can  be  varied,  depending  on  the  length 
of  time  following  installation.  Excess  pore  pressures  and  reduced  shear 
strengths  exist  for  periods  up  to  30  days  or  longer. 

Total  Stress  Methods 


a-method 

11.  The  a-method  (API  1978,  Vesic  1972)  is  based  on  observations  that 
the  ratio  of  the  pile-soil  adhesion,  ca  ,  to  the  undrained  shear  strength  of 
the  soil  decreases  as  the  undrained  shear  strength  increases.  The  ratio  is  a 
parameter  a  that  is  given  as  a  function  of  the  undrained  shear  strength. 

The  total  pile  capacity,  Qu  ,  is  then  computed  by  using  Equation  1: 


where 


su  =  undrained  shear  strength  along  the  side  of  the  pile  or  at  the  pile 
tip 

As  =  area  of  the  pile  surface  equal  to  the  perimeter  of  the  pile  multi¬ 
plied  by  the  pile  length 

Nq  =  cohesive  bearing  factor  (for  piles  in  clay,  Nc  can  be  taken  as  9 1 

At  =  cross-sectional  area  of  the  pile  tip 

In  practice,  values  of  a  have  been  computed  using  the  averages  of  ca  and 
su  for  the  total  embedded  lengths,  and  the  side-frictional  capacity  is  com¬ 
puted  by  dividing  the  subsurface  into  a  number  of  layers  using  the  average  un- 
drained  shear  strength  of  each  layer. 

12.  There  are  factors  that  influence  the  accuracy  of  the  a-method 
Pile-soil  adhesion  is  not  only  dependent  on  the  undrained  shear  strength  but 
is  also  dependent  on  the  soil  composition,  stress  history,  and  sensitivity. 

It  has  also  been  shown  (Kraft,  Focht,  and  Amarasinghe  1981)  that  there  is  also 
some  dependence  of  adhesion  on  the  depth  of  embedment  and  the  soil  stiffness 
rather  than  being  limited  to  values  prescribed  by  a  . 

13.  In  spite  of  these  shortcomings,  the  a-method  has  been  widely  used 
and  will  continue  to  be  used.  More  reliable  results  would  be  obtained  if  un¬ 
published  load-test  results  from  short,  onshore  piles  were  analyzed.  Other 
investigators  (Dennis  and  Olson  1983)  (Kraft,  Focht,  and  Amarasinghe  1981) 
have  conducted  such  studies  for  steel  pipe  piles  for  offshore  purposes. 

14.  The  literature  contains  a  number  of  recommendations  to  obtain  pile- 
soil  adhesion,  according  to  the  a-method.  Figure  1  from  Tomlinson  (1957) 
shows  curves  for  adhesion  as  a  function  of  shear  strength.  Table  1  contains 
values  of  adhesion  recommended  by  Tomlinson  (1964).  Figure  2  from  Vesic 
(1977)  illustrates  values  of  adhesion  from  a  number  of  tests  including  tests 
on  cast-in-place  concrete  piles.  Figure  3  (Tomlinson  1971)  gives  adhesion 
factors  for  stiff  clays  for  piles  with  a  penetration  ratio,  depth  of  penetra¬ 
tion  in  clay  divided  by  pile  diameter,  greater  than  20  for  piles  other  than 
H-piles.  A  recommendation  is  to  use  a  =  0.4  for  piles  with  a  penetration 
ratio  between  8  and  20.  Figure  4  indicates  recommendations  by  other  re¬ 
searchers  (Vesic  1977). 

15.  Results  of  statistical  analyses  by  other  researchers  have  been  pub¬ 
lished  recently.  Kraft,  Focht,  and  Amarasinghe  (1981)  analyzed  the  results  of 
testing  on  closed  and  open-end  pipe  piles,  timber,  and  concrete  piles.  Values 
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Figure  2.  Comparison  of  side  resistance  and  undrained 
shear  strength  (Vesic  1977) 
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Undrained  shear  strength,  (c  ),  lb/ft 
Figure  3-  Adhesion  factors  for  stiff  clays  (Tomlinson  1970 


of  a  were  computed  by  a  method  proposed  by  the  API  (1978)  and  outlined  as 
follows.  For  highly  plastic  clays  with  liquid  limits  and  plasticity  indices 
greater  than  50  and  35  percent,  respectively,  a  =  1  except  for  overconsoli¬ 
dated  clays  where  c_  must  be  less  than  1  ksf  or  the  undrained  shear  strength 

of  the  soil  if  it  was  normally  consolidated  clay,  (su)  .  For  other  clays, 

nc 

a  varies  linearly  from  1  at  su  =  0.5  to  0.5  at  su  =  1.5  ksf. 

16.  Computed  values  of  a  were  plotted  as  a  function  of  pile  length, 

L  ,  and  a  parameter  relating  pile  stiffness  to  soil  stiffness. 
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D  =  pile  diameter 
*max  =  nax^nuja  side  resistance 

A  =  cross-sectional  area  of  the  pile 


(2) 


10 


c*  »  kN/m 


\ 


3 1  <0 
o  |  o 


Undrained  cohesion,  (c^),  kips/ft 

Figure  4.  Adhesion  factors  for  piles  in  clay  (Vesic  1977) 

E  =  Young's  modulus  of  the  pile 
u*  =  vertical  pile  movement  at  maximum  side  resistance 
The  relationships  for  o  for  normally  consolidated  soils  are 

a  =  1.486  -  0.126  ln(L)  (3) 

a  =  1.024  -  0.070  ln(*3)  (4) 

Similar  expressions  were  obtained  for  overconsolidated  soils.  The  linear 
correlation  coefficients  are  -0.27  and  -0.33  for  Equations  3  and  4,  respec¬ 
tively.  The  authors  also  determined  the  ratios  of  computed  to  measured  capac¬ 
ity.  A  mean  of  1.10  and  standard  deviation  of  0.34  was  obtained  by  using 
Equation  3(  and  a  mean  of  1.09  and  standard  deviation  of  0.31  was  obtained  by 
using  Equation  4.  So,  in  spite  of  the  fact  that  the  expressions  do  not  corre¬ 
late  very  well,  reasonably  accurate  predictions  of  capacity  can  be  made. 
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17.  Predictions  of  capacity  were  obtained  by  using  three  different 
{■•methods  for  steel  pipe  piles  in  normally  or  lightly  overconsolidated  clay 
(Gardner  1977).  The  first  method  used  o  as  proposed  by  Tomlinson  (1957)  for 
steel  piles.  A  second  method  for  layered  soil  systems  used  values  of  a  pro¬ 
posed  by  Tomlinson  (1971),  and  the  third  is  the  method  proposed  by  API  (1978). 
Ratios  of  the  computed  to  the  measured  capacities  for  57  compression  tests 
were  calculated.  The  respective  mean  and  standard  deviations  for  the  ratios 
were  0.74  and  0.48,  1.13  and  0.45,  1.07  and  0.37  for  each  of  the  three  methods. 
A  new  method  was  proposed  as  follows: 


Q  =  ac  F  F.  A  +  9c  F  A. 
u  ucLs  uct 


(5) 


where 

cu=  average  su 

Fc  =  ratio  of  the  undrained  shear  strength  using  unconsolidated  un¬ 
drained  triaxial  tests  to  the  strengths  obtained  by  other  methods 

FL  =  1.0  for  lengths  up  to  100  ft  and  varies  linearly  to  1.8  at  175  ft 

a  =  1.0  for  values  of  c  F  of  600  psf  or  less  and  varies  linearly  to 

u  c  _ 

0.5  for  c  F  equal  to  1,200  psf  and  to  0.3  for  c  F  equal  to 
u  c  u  c 

5,000  psf 

The  mean  for  the  ratios  of  computed  to  measured  capacities  was  1.00  and  the 
standard  deviation  was  0.22. 

18.  In  view  of  the  previous  discussion,  the  recommendation  is  to  deter¬ 
mine  a  according  to  the  method  proposed  by  Dennis  and  Olson  (1983)  for  piles 

100  ft  or  less.  As  shown  in  Figure  5,  the  recommendation  for  a  is  within 

the  range  of  values  of  a  that  have  been  recommended  by  others.  It  is  also 

recommended  that  the  undrained  shear  strength  be  correlated  to  the 

unconsol idated-undrained  triaxial  strength.  It  is  noted  here  that  the  use  of 
Equations  3  and  4  could  lead  to  overpredictions  of  capacity  for  piles  less 
than  100  ft. 
p-method 

19.  The  p-method  predicts  the  side  resistance,  f  ,  as  a  function  of 
the  in  situ  vertical  overburden  pressure,  o' 


2 

Undrained  cohesion,  (cu>,  kips/ft 

Figure  5.  Adhesion  factors  as  recommended  by  Dennis  and  Olson  (1983) 

pressure  and  the  tangent  of  the  effective  angle  of  internal  friction,  *’  ,  as 
shown  in  Equation  7. 

6  =  (1  -  sin  *')  tan  (0CR)1/2  (7) 

where  OCR  is  the  overconsolidation  ratio.  Another  recommendation  is  to  de¬ 
termine  b  as  a  function  of  pile  length  for  normally  and  lightly  overconsoli¬ 
dated  clays.  Figure  6  (Meyerhof  1976)  indicates  that  the  trend  is  for  8  to 
decrease  with  increasing  pile  length,  particularly  for  piles  longer  than 
100  ft.  Meyerhof  (1976)  also  recommends  that  s  for  overconsolidated  soils 
be  taken  as 

i  /? 

8  =  1.5  (1  -  sin  ♦')  tan  (OCR)  (8) 

20.  Statistical  analyses  of  the  8-method  were  also  presented  by  Kraft, 


Depth  of  penetration,  ft 


The  linear  correlation  coefficient  for  Equations  9  and  10  were  -0.58  and 
-0.71,  respectively.  The  respective  mean  and  standard  deviation  of  the  ra¬ 
tio  of  computed  to  measured  capacity  were  1.05  and  0.22  for  Equation  9,  and 
1.03  and  0.19  for  Equation  10.  Values  of  8  did  not  correlate  for  overcon¬ 
solidated  soils.  A  recommendation  for  the  use  of  the  8-method  is  given 
subsequently . 

1-method 

21.  The  x-method  correlates  the  theoretical  passive  earth  pressure  to 
the  side  resistance  (Vijayvergiya  and  Focht  1972). 

f  =  l( o'  ♦  2c  )  (11) 

m  m 

where 

o'  =  mean,  effective,  overburden  pressure 

ID 

c,,,  =  mean,  undrained  shear  strength 

both  for  the  length  of  embedment.  Figure  7  is  a  plot  of  X  versus  pile 
length  for  pipe  piles.  The  figure  indicates  that  there  is  a  strong 
correlation  between  the  factor  X  and  pile  embedment,  that  it  decreases  with 
length . 

22.  Results  similar  to  the  analyses  of  the  previous  two  methods  were 
also  obtained  for  the  x-method.  Kraft,  Focht,  and  Amarasinghe  (1981)  deter¬ 
mined  relationships  for  X  as  a  function  of  pile  length  and  *3  for  normally 
consolidated  soils. 


X  =  0.296  -  0.032  ln(L) 

(12) 

X  =  0.178  -  0.016  ln(*3) 

(13) 

23.  The  expressions  for  overconsolidated  soils,  soils  having  ratios  of 
the  mean,  undrained  shear  strength  to  the  mean,  effective,  overburden  pressure 
greater  than  0.4,  are 


x  =  0.488 

-  0.078  ln(L) 

(14) 

X  =  0.23 2 

-  0.032  ln(w3) 

(15) 

0 


i 

I 

factor  necessary  is  (Kraft,  Focht,  and  Amarasinghe  1981)  where 

(16) 

v° 

For  the  a-method  and  the  X -method,  respectively,  8  then  would  be 

i 

i 

8  =  a81  (17)  | 

! 

8  =  (1  +  28,)  (18)  | 

I  I 

( 

Figure  8  shows  a  comparison  of  the  three  methods  using  8^  =  0.2  .  The  re¬ 
sults  obtained  for  the  a-method  and  the  X-methods  would  result  in  lower  pre¬ 
dicted  capacities  in  this  particular  case,  but  not  for  cases  when  8^  >  0.2  .  j 

The  curves  from  the  8-method  represent  mean  values  for  soils  for  all  values  of  ! 

Effective  Stress  Methods 


26.  The  development  of  various  effective  stress  methods  as  a  means  of 
calculating  pile  capacity  has  occurred  as  a  result  of  researchers’  different 
attempts  to  calculate  the  state  of  stress  in  the  soil  around  piles  during 
loading.  The  state  of  stress  is  updated  to  represent  conditions  in  the  soil 
after  pile  installation,  after  soil  consolidation,  and  at  pile  failure.  Solu¬ 
tions  require  theoretical  models  to  represent  the  strain  behavior  of  soil  that 
occurs  during  each  of  the  conditions  mentioned  above,  and  a  soil  model  to  rep¬ 
resent  stress-strain-strength  relationships.  The  accuracy  of  any  effective 
stress  method  is,  thus,  dependent  on  the  ability  to  model  the  strains  that  oc¬ 
cur  and  the  stress-strain  behavior  of  the  soil. 

Modeling  pile  installation 

27.  Pile  installation  has  been  modeled  using  cylindrical  cavity  expan¬ 
sion  (Esrig  et  al.  1977;  Esrig  and  Kirby  1979;  Heydinger  and  O'Neil,  in  prepa¬ 
ration;  Kirby,  Esrig,  and  Murphy  1983;  Kirby  and  Roussel  1979  and  1980;  Kirby 
and  Wroth  1977),  or  other  equivalent  methods  (Carter,  Randolph  and  Wroth  1979; 
Randolph,  Carter,  and  Wroth  1979;  Wroth,  Carter,  and  Randolph,  1979).  Accord¬ 
ing  to  the  methods,  it  is  assumed  that  the  soil  only  moves  outward,  radially 
(direction  perpendicular  to  the  pile  surface),  resulting  in  plane  strain 
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Figure  8.  Comparison  of  predictive  methods 

conditions.  For  cylindrical  cavity  theory,  the  soil  is  assumed  to  behave  as 
an  elastic,  perfectly  plastic  material.  For  the  other  formulations,  work¬ 
hardening  or  softening  effects  can  occur.  Concepts  of  critical  state  soil 
mechanics  (Kirby  and  Wroth  1977,  Schofield  and  Wroth  1968)  for  soils  at  large 
strains  are  included  in  the  models.  Accordingly,  soil  reaches  a  unique  rela¬ 
tionship  of  stress-strain-volume  that  depends  on  the  initial  state  of  stress 
and  strain  and  the  stress  history.  The  theories  are  used  to  compute  the  state 
of  stress  in  the  soil  and  the  excess  pore  pressure  distribution  immediately 
following  pile  installation. 


Modeling  soil  consolidation 

28.  Existing  consolidation  models  assume  that  pore  water  movements  oc¬ 
cur  only  in  the  radial  direction.  Initially,  it  was  assumed  that  the  soil 
would  remain  at  the  critical  state  of  effective  stress  during  consolidation 
(Esrig  et  al.  1977,  Kirby  and  Esrig  1979).  Other  solutions  used  an  elastic 
soil  model  (Randolph  and  Wroth  1979),  an  elastic  soil  model  with  stiffness 
that  increased  with  distance  from  the  pile  surface  (Heydinger  and  O'Neill,  in 
preparation;  Lelfer,  Kirby,  and  Esrig  1979),  an  elasto-plastic  soil  model 
(Kavvadas  and  Baligh  1982;  Miller,  Murff,  and  Kraft  1978)  and  rigorous  formu¬ 
lations  that  modeled  the  soil  as  an  elasto-plastic  material  (Carter,  Randolph, 
and  Wroth  1979;  Randolph,  Carter,  and  Wroth  1979;  Wroth,  Carter,  and  Randolph 
1979).  Published  results  for  the  soil  adjacent  to  the  pile  give  the  changes 
in  mean  effective  stress  or  radial  effective  stress  as  a  function  of  the  ex¬ 
cess  pore  pressure  immediately  after  installation.  It  is  generally  accepted 
(Kirby,  Esrig,  and  Murphy  1983;  Randolph,  Carter,  and  Wroth  1979;  Wroth,  Car¬ 
ter,  and  Randolph  1979)  that  soil  will  return  to  KQ  ,  conditions  for  normally 
consolidated  soil,  with  the  radial  effective  stress  equal  to  the  maximum  prin¬ 
cipal  stress. 

Pile  loading 

29.  Various  attempts  have  been  made  to  determine  the  effects  of  pile 
loading  on  the  state  of  stress  in  the  soil.  Initially,  it  was  assumed  that 
the  soil  was  at  the  critical  state  after  consolidation  so  that  the  mean  effec¬ 
tive  stress  would  remain  constant  during  shearing  (Esrig  et  al.  1977).  Other 
investigators  (Kraft  1982),  relying  on  findings  that  the  soil  would  not  be  at 
the  critical  state,  have  used  stress  paths  to  predict  the  state  of  stress  at 
pile  failure.  Pile  loading  effects  have  also  been  investigated  by  using  the 
finite  element  method  with  a  one-  or  two-dimensional  scheme  (Baguelin  and 
Frank  1980;  Baguelin,  Frank,  and  Jezequel  1982;  Potts  and  Martins  1982).  Pile 
loading  has  also  been  modeled  with  a  three-dimensional  finite  element  program 
using  the  state  of  stress  after  consolidation  as  input,  but  which  did  not 
model  the  effective  stress  changes  during  pile  loading  (Heydinger  and  O’Neill, 
in  preparation). 

Pile  capacity  predictions 

30.  Effective  stress  methods  predicting  side  resistance  have  been  pro¬ 
posed.  Esrig  and  Kirby  (1979)  developed  a  set  of  curves  to  determine  6  as 
a  function  of  the  soil  OCR  and  plasticity  index.  These  curves  result  in 


overpredictions  of  capacity,  with  s  ranging  from  0.3  from  a  normally  con¬ 
solidated  soil  to  2.0  for  highly  overconsolidated  soils. 

31.  A  large  research  effort  funded  by  a  number  of  organizations  has  re¬ 
sulted  in  at  least  four  effective  stress  methods  (Kraft  1982).  The  four  meth¬ 
ods  reported  by  Kraft  utilize  different  approaches  to  predict  capacity.  Ac¬ 
cording  to  the  first  two  methods  reported  by  Kraft  (1982),  the  soil  is  at  the 
critical  state  after  installation,  but  is  not  for  the  other  two  methods.  The 
ratio  of  the  change  in  mean  effective  stress  during  consolidation  to  the  ex¬ 
cess  pore  pressure  varies  from  0.25  to  about  0.8,  depending  on  the  method. 

The  first  and  third  method  listed  by  Kraft  (1982)  assumed  that  the  mean  effec¬ 
tive  stress  remains  constant  during  pile  loading.  However,  it  is  assumec  that 
the  mean  effective  stress  decreases  for  the  other  methods.  The  comparisons  of 
predicted  and  measured  capacities,  using  the  four  effective  stress  methods  and 
the  a  and  X  methods,  indicated  that  the  effective  stress  methods  could  be 
used  as  accurately  as  the  a  and  X  methods.  An  underlying  problem  with  the 
effective  stress  methods  is  in  determining  soil  parameters  that  are  used  to 
compute  the  changes  in  stress. 

32.  Two  other  effective  stress  methods  based  on  predictions  of  the 
effective  radial  stress  after  consolidation  have  been  proposed.  Results  from 
a  three-dimensional  finite  element  program  indicate  that  the  total  radial 
stress  acting  on  the  pile  does  not  change  significantly  during  pile  loading 
(Heydinger  and  O'Neill,  in  preparation).  Therefore,  by  neglecting  the 
hydrostatic  pore  pressure,  the  effective  radial  stress  after  consolidation  can 
be  used  for  the  total  radial  stress  in  the  following  expression: 


f  =  a  tan  $  „ 
rc  ss 


where  ♦  is  the  angle  of  friction  between  the  pile  and  the  soil.  The  rec- 

SS 


ooanended  equation  for  side  shear  for  the  one  method  is 


( 1  +  sin  $ ')  Ac' 

_  cs  rc  , 

f  =  - : - : -  s  +  -  u  tan  4> 

sin  <>'  u  u  e  ss 

cs  cs  e 


where 


=  effective  angle  of  friction  for  soil  at  the  critical  state 


(large  shearing  strains) 

=  undrained  shear  strength  for  soil  at  the  critical  state 


( Ao *  /u  )  =  change  in  effective  radial  stress  divided  by  the  excess 
rc  e  pore  pressure,  u 


The  equation  for  u  is 


/r  \1/2 

k  -  *  2(su)cs 


(21) 


where 


G  =  undrained  shear  modulus 


P' 

*cs 
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33.  The  other  method  used  the  results  from  two  consolidation  models 
(Carter,  Randolph,  and  Wroth  1979;  Heydinger  and  O'Neill,  in  preparation)  to 
obtain  the  effective  radial  stress  after  consolidation.  An  expression  for  the 
stress  was  obtained  from  the  results  of  a  parametric  study. 


—  =  4.80  -  4.57  log  M 
s 

u 


(22) 


where  M  =  (6  sin  <*')/(3  -  sin  $') 
taken  as  the  in  situ  value. 

34.  Estimates  of  8( 


The  undrained  shear  strength  should  be 


f/4»^Q)  can  be  obtained  from  Equations  20  and  22. 


A  range  of  values  of  <>'  , 

CS 


(M 


/o' 


cs 


VO 


( Ao  /u  )  and  G/s,,  were  substi- 


rc 


tuted  into  Equations  20  and  21,  resulting  in  a  range  of  values  of  8  between 
0.3  and  0.6.  A  similar  procedure  was  followed  with  Equation  22  resulting  in  a 
range  of  values  of  8  between  0.3  and  0.55.  It  is  apparent  that  the  computed 
capacity  can  be  overestimated  unless  accurate  estimates  of  the  respective  pa¬ 
rameters  are  made. 


Comparison  of  Predictions 

35.  Four  instrumented  pile-load  tests  were  used  to  compare  predic¬ 
tions  for  side  shear.  The  four  pile-load  tests  are  from  insensitive,  over¬ 
consolidated  clay  at  the  University  of  Houston  Central  Campus  (UHCC)  (O'Neill, 
Hawkins,  and  Mahar  1 98 1 ) ;  slightly  sensitive,  normally  consolidated  clay  in 
San  Francisco  Bay  Mud  (SFBM)  (Kirby  and  Roussel  1979  and  1980);  sensitive, 
lightly  overconso L idated  clay  at  St.  Alban  near  Quebec  (Bjerrum  and  Simons 
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I960;  Konrad  1977;  Roy  et  al.  1981);  and  insensitive,  normally  consolidated 
clay  at  the  East  Atchafalya  Basin  Protection  Levee  (EABPL)  (USAE  District,  New 
Orleans  1977;  Shilstone  Testing  Laboratory  1976).  Mea3ured-side-shear  and 
computed-side-shear  distributions  are  tabulated  and  ratios  of  the  computed  to 
the  measured  side  shear  are  given. 

36.  The  side  shear  was  computed  at  different  depths  for  each  of  the 
four  tests  using  five  methods.  For  the  a-method,  Equation  5  was  used  as  rec¬ 
ommended  by  Dennis  and  Olson  (1983).  Side  shear  was  computed  for  the  B-method 
according  to  Equation  8,  obtained  from  Meyerhof  (1976).  Equations  12  and  14 
from  the  statistical  studies  of  Kraft,  Focht,  and  Amarasinghe  (1981)  were  used 
to  compute  x  for  the  X-method.  The  other  two  methods  were  obtained  using 
Equations  20  and  22.  The  parameters  for  Equation  20  were  estimated  for  the 
St.  Alban  and  EABPL  tests  since  they  are  not  available. 


UHCC  test 


37.  The  test  at  the  UHCC  site  was  part  of  a  testing  program  sponsored 
by  the  Federal  Highway  Administration.  The  stratigraphy  for  the  UHCC  test 
consists  of  1.5  ft  of  clay  fill  underlain  by  two  formations  of  clay  which  were 
preconsolidated  by  dessication.  The  top  4  ft  of  the  lower  formation,  from  a 
depth  of  26  to  30  ft,  was  appreciably  softer  than  the  soils  above  and  below. 
Undrained  shear  strengths  were  obtained  from  unconsolidated,  undrained  (UU) 
and  consolidated,  isotropic,  undrained  (CIU)  triaxial  testing.  The  test  pile 
and  test  that  was  modeled  for  this  research  was  reference  Pile  1  and  Test  1 , 
conducted  19  days  after  pile  driving.  The  test  pile  was  a  steel  pile, 
10.75-in.  outside  diameter  (0D)  and  10-in.  inside  diameter  (ID).  The  pile  was 
driven  43  ft  below  ground  surface  and  was  instrumented  with  strain  gage  loca¬ 
tions  5  ft  apart.  The  pile  was  loaded  in  increments  of  approximately  10  tons 
that  were  maintained  for  periods  of  about  1  hr  until  plunging  failure  occur¬ 
red.  The  measured-  and  computed-side  shear  for  the  test  i3  shown  in  Table  2. 


SFBM  test 


38.  The  SFBM  test  was  conducted  at  Hamilton  Air  Force  Base.  The  sub¬ 
surface  conditions  consist  of  a  single  deposit  of  relatively  homogeneous  ma¬ 
rine  clay  to  a  depth  of  at  least  50  ft.  It  is  a  very  soft  clay  except  for  the 
top  6  ft  which  is  dessicated.  The  undrained  shear  strengths  were  obtained 
from  unconfined  compression,  UU  and  field-vane  shear  tests.  The  test  pile  was 
4.5-in.  0D  by  4.0-in.  ID  sections  Jacked  to  a  depth  of  40  ft.  Instrumented 
with  strain  gage  locations  spaced  at  5-ft  intervals,  the  pile  was  loaded  to 


EMWW 


»*<  *  - 


M* 

m 


w 


V  V 
•>SV 


55 


Xy^  v 

V  V  V 


failure  in  approximately  10  min  by  applying  30  small-load  increments  every 
20  sec.  The  measured  and  computed  side  shears  are  shown  in  Table  3. 

St.  Alban  test 

39.  The  St.  Alban  test  was  conducted  near  Quebec,  Canada.  The  subsur¬ 
face  consists  of  a  foot  of  topsoil,  4  ft  of  weathered  clay  crust,  27  ft  of 
soft  silty  clay  of  marine  origin,  13  ft  of.  very  soft  clayey  silt,  and  a  dense 
layer  of  sand  below  a  depth  of  45  ft.  The  soil  was  overconsolidated  due  to 
secondary  consolidation  effects.  The  undrained  shear  strengths  were  obtained 
from  the  average  of  field  vane  tests  which  were  less  than  strengths  obtained 
from  UU  tests.  For  this  report,  Pile  6  was  used  to  compute  side  resistance. 
The  pile  was  a  steel  casing  with  an  8.625-in.  0D  and  a  7.99-in.  ID.  The  pile 
was  Jacked  to  a  depth  of  25.3  ft  below  ground  surface  with  an  oversized  pilot 
hole  4  ft  deep.  Loads  were  measured  only  at  the  top  and  the  bottom  of  the 
pile,  leaving  the  average  side  friction  for  use  in  the  comparisons.  Results  of 
a  test  conducted  at  472  hr  were  used.  Load  increments  of  about  139  lb  were 
maintained  for  periods  of  about  15  min  until  failure  occurred  after  3  hr.  The 
measured  and  computed  side  shears  are  given  in  Table  4. 

EABPL  test 

40.  The  load  test  was  on  the  EABPL,  Test  1-2.  The  soil  profile  con¬ 
sists  of  84  ft  of  soft  to  stiff  clay  underlain  by  20  ft  of  silt.  Undrained 
shear  strengths  were  estimated  using  unconfined  compression,  UU  and  CIU  tests. 
A  concrete  prestressed  pile  14  sq  in.  was  used.  Instrumented  gage  locations 
were  placed  at  8-ft  intervals.  The  pile  was  driven  and  then  tested  approxi¬ 
mately  one  month  after  driving.  The  pile  was  loaded  in  increments  of  10  tons 
for  periods  of  about  1  hr  except  for  the  40-ton  load  which  was  held  for  24  hr. 
The  comparisons  of  side  shear  are  given  in  Table  5. 

Conclusions  and  Recommendations 

41.  Predictions  of  axial  capacity  can  be  obtained  most  reliably  by  us¬ 
ing  total  stress  methods.  The  recommendation  is  to  use  the  lesser  of  the  ca¬ 
pacities  computed  with  the  a-method  by  Dennis  and  Olsen  (1983)  shown  in  Fig¬ 
ure  5,  or  the  0-method  using  Equation  8.  The  predicted  side  capacity,  as  it 
is  developed  along  piles,  can  vary  significantly  for  the  two  methods.  The 
computed  side  shear  is  greater  near  the  ground  surface  for  the  a-method  and 
lesser  at  greater  depths,  compared  to  the  e-method.  The  e-method  is 


consistent  with  observations  that  the  unit  side  shear  increases  with  depth. 

The  X-method  can  result  in  unconservative  values  of  capacity.  The  analyses 
for  the  recommended  procedures  were  determined  principally  from  the  result  of 
pile  pipes  so  additional  data  with  concrete  piles  and  H-piles  are  desirable. 

42.  Effective  stress  methods  generally  result  in  unconservative  predic¬ 
tions  of  pile  capacity.  There  are  too  many  unknown  soil  parameters  that  are 
necessary  for  effective  stress  methods.  Additionally,  the  methods  were  devel¬ 
oped  for  impermeable  cylindrical  piles  and  not  for  concrete  or  H-piles.  Other 
effects  of  driving  that  the  models  do  not  account  for  are  lateral  pile  wobble 
and  repeated  soil  shearing  during  pixe  driving. 


PART  III:  PILE  SIDE  SHEAR-DISPLACEMENT  FUNCTION 


Introduction 


43.  Research  on  pile  side  shear-displacement  functions  has  occurred  as 
a  result  of  experimental  observations  that  side  shear  that  is  mobilized  is  de¬ 
pendent  on  vertical  pile  displacement  (Coyle  and  Reese  1966,  D’Appolonia, 
Poulos,  and  Ladd  1971)  and  because  of  their  use  in  soil-structure  interaction 
programs  (Radhakrishnan  and  Parker  1975).  Piles  have  been  instrumented  with 
strain  gages  at  various  levels  in  order  to  measure  the  distributions  of  loads 
along  piles.  Load  distribution  curves  were  then  used  to  compute  mobilized 
side-shear  and  vertical  pile-displacement  relationships.  The  relationships, 
f-z  curves,  are  nonlinear  functions  of  side  shear  and  vertical  pile  displace¬ 
ment.  Soil-structure  interaction  programs  are  algorithms  that  model  piles  as 
discrete  elements.  The  programs  use  f-z  curves  to  define  the  load-transfer 
behavior  at  the  pile-soil  interface  in  order  to  determine  load-deformation  be¬ 
havior  of  piles.  The  purpose  of  this  section  is  to  discuss  recommendations 
for  the  determining  f-z  curves. 

Methods  for  Computing  f-z  Curves 

44.  Four  methods  for  computing  f-z  curves  were  investigated.  The 
computational  procedures  are  presented  in  this  section.  The  next  section  con¬ 
sists  of  an  evaluation  of  the  four  methods. 

Coyle  and  Reese  method 

45.  Coyle  and  Reese  (1966)  used  the  results  of  full-scale  and  labora¬ 
tory  testing  in  order  to  develop  a  procedure  for  determining  f-z  curves. 

The  recommendation  was  to  use  Figure  9  to  determine  the  effective  shear 
strength  and  Figure  10  to  determine  the  load  transfer,  normalized  by  the  ef¬ 
fective  shear  strength,  for  a  number  of  pile  movements.  The  load  transfer  is 
computed  by  multiplying  the  effective  shear  strength  by  the  normalized  load 
transfer  for  a  number  of  vertical  pile  movements. 

Vljayverglya  method 

46.  A  parabolic  expression  was  proposed  by  Vijayvergiya  (1977b)  to  com¬ 
pute  f-z  curves.  According  to  this  method,  the  pile  movement  should  be  nor¬ 
malized  by  the  critical  pile  movement,  zc  ,  that  occurs  when  the  maximum  side 
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Figure  9.  Shear  strength  reduction  factor,  K 
(Radhakrishnan  and  Parker  1975) 
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,  is  mobilized.  The  equation  is 
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The  recommendation  was  to  use  z_  equal  to  0.25  in.  fm_„ 
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is  the  side  shear 


capacity  which,  for  practical  purposes,  can  be  determined  by  any  method. 


Kraft,  Ray,  and  Kagawa  method 


47.  The  method  proposed  by  Kraft,  Ray,  and  Kagawa  (1981)  is  a  modifica¬ 
tion  of  an  approximate  elastic  solution  proposed  by  Randolph  and  Wroth  (1978). 
According  to  the  method,  the  vertical  displacements  at  a  given  level  along  the 
pile  are  approximated  by  a  curved  surface,  with  displacements  decreasing  as 
the  radial  distance  from  the  center  of  the  pile  increases.  The  curved  surface 
is  approximated  with  concentric  cylinders  of  shear.  It  is  assumed  that  the 
radial  displacements  are  negligible  compared  to  the  vertical  displacement, 
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Figure  10.  Normalized  load  transfer  versus  vertical  pile 
movement  (Radhakrishnan  and  Parker  1975) 

which  then  results  in  simple  shear  conditions.  For  simple  shear,  the  product 
of  the  shear  stress,  t  ,  at  radial  distance,  r  ,  from  the  pile  center  and  the 
radial  distance  is  constant.  It  is  assumed  that  the  shear  stresses  are  negli¬ 
gible  at  a  radial  distance,  rm  ,  where  the  boundary  of  the  zone  of  influence 
is  given  in  Equation  24, 

rffl  =  2.5  lo  (1  -  v)  (24) 

where 

L  =  pile  length 
v  =  Poisson’s  ratio  of  the  soil 

p  =  ratio  of  the  shear  modulus  at  depth  L/2  to  the  shear  modulus  at 
the  pile  tip. 
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48.  The  load -displacement  relationship  is  obtained  by  integrating  the 
shear  strains  from  the  pile  surface,  r0  ,  to  the  boundary  of  the  zone  of  the 
Influence.  Since  the  product  r  times  t  is  constant,  Equation  25  is  ob¬ 
tained  for  the  relationship, 
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dr 
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where 

zfl  s  vertical  displacement  at  the  pile  surface 
s  soil  shear  stress  at  the  pile  surface 
G  =  soil  shear  modulus 

Because  of  the  effects  of  pile  installation  and  soil  console  ion,  the  shear 
modulus  increases  approximately  linearly  with  radial  distance  from  the  pile 
until  the  undrained  modulus  is  reached  at  a  radius  of  rm  .  It  was  shown  by 
Kraft,  Ray,  and  Kagawa  (1981)  that  the  average  shear  modulus  after  soil  con¬ 
solidation  does  not  vary  significantly  from  the  initial  shear  modulus  for  un¬ 
disturbed  soil,  ,  for  soil  at  very  low  strains.  Therefore,  the  radial 
variation  of  the  shear  modulus  due  to  pile  installation  was  neglected  and  the 
variation  of  the  shear  modulus  due  to  shear  stress  was  considered. 

49.  A  secant  modulus  formulation  was  used  to  approximate  the  variation 
of  the  shear  modulus  with  shear  stress.  The  variation  is  expressed  by  the 
following  hyperbolic  equation 


Go  =  Gi 


where 

Rf  =  stress-strain  curve-fitting  constant 
Tmax  =  maximuin  shear  stress  that  is  mobilized  at  pile  failure 
No  specific  recommendations  were  given  for  estimating  R^  when  shear  stress- 
strain  data  are  not  available.  It  is  typically  in  the  range  of  0.9  to  1.0. 

50.  An  expression  for  the  vertical  pile  displacement  is  obtained  by 
substituting  Equation  26  for  Equation  25.  The  resulting  expression  is 
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where 


Tmax 


(28) 


Equation  27  can  be  solved  for  a  number  of  vertical  pile  displacements  by  se¬ 
lecting  different  values  of  shear  stress  that  are  less  than  the  maximum  shear 
stress.  The  resulting  data  can  be  plotted  as  a  curve  of  shear  stress  versus 
vertical  pile  movement,  which  would  be  the  desired  f-z  curve. 

Heydinger  and  O'Neill  method 

51.  A  general  equation  for  side  shear  as  a  function  of  pile  movement 
normalized  by  pile  diameter,  z/D  ,  was  proposed  by  Heydinger  and  O'Neill  (in 
preparation).  The  equation  for  side  shear  is 


Efz  =  slope  of  the  initial  tangent  to  a  curve  of  f  versus  z/D 
m  =  shape  parameter 

For  the  investigation,  three  well- instrumented  pile-load  tests  were  modeled 
extensively  using  an  approximate  solution  to  determine  the  effects  of  pile  in¬ 
stallation  and  a  finite  element  formulation  on  model  pile  loading.  A  parame¬ 
tric  study  was  conducted  with  the  finite  element  program  in  order  to  determine 
the  effects  of  pile  properties  on  the  shapes  of  f-z  curves.  The  computed  f-z 
curves  were  optimized  for  the  parameters  in  terms  of  soil  and  pile  properties. 

52.  Expressions  for  the  parameters  E^z  and  m  were  obtained  as  fol¬ 
lows.  Efz  is  given  as  the  ratio  of  the  initial  undrained  soil  modulus,  Eu  , 
and  a  parameter  K  .  An  expression  for  K  was  determined  as  a  function  of 
the  pile  length  to  diameter  ratio,  L/D  as  shown  in  Figure  11.  The  expres¬ 
sion  is 
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The  shape  parameter  was  plotted  as  a  function  of  the  pile  L/D  ratios  as 
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Figure  11.  In  K  versus  In  (L/D) 


shown  in  Figure  12.  The  intercept  values  for  m  ,  referred  to  as  mQ  ,  were 
plotted  as  a  function  of  the  ratio  / Eu\  jpa  as  shown  in  Figure  13,  where 


(Eu)av/S 


is  the  average  undrained  soil  modulus  for  a  soil  stratum  and  Pa  is 


atmospheric  pressure.  A  linear  regression  analysis  was  used  to  obtain  Equa¬ 
tion  31  for  the  shape  parameter. 


m  =  exp  <0.12  +  0.54  In 
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The  reason  that  the  pile  movements  were  normalized  by  the  pile  diameters  and 


(E)  was  normalized  by  atmospheric  pressure  was  to  nondimensional ize  the 
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Figure  12.  In  m  versus  In  (L/D) 

Example  of  Predictive  Methods 

53.  The  following  computational  results  illustrate  the  use  of  the  four 
predictive  methods.  For  this  example,  a  pile  with  an  outside  diameter  of  1  ft 
and  a  length  of  75  ft  was  modeled.  The  f-z  data  was  computed  for  a  depth  of 
22.5  ft. 

54.  Soil  properties  representative  of  a  normally  consolidated  soil  were 
used  for  the  computations.  A  wet  unit  weight  of  112.4  pcf  and  a  depth  to  the 
water  table  of  15  ft  were  selected.  The  mean  effective  overburden  pressure 
was  computed  assuming  kQ  equal  to  0.58.  An  undrained  shear  strength  of 

721  psf  was  computed  using  the  ratio  3u/(3y0  equal  to  0.35  (su/p^  =  0.49)  . 
The  maximum  side  shear  was  computed  using  the  a-method  as  recommended  by 
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Figure  13.  In  m  versus  In 
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Dennis  and  Oison  (1983).  Thus  a  =  0.90  and  fmax  =  648  psf  were  obtained. 
The  undrained  soil  modulus,  equal  to  865,800  psf,  was  computed  using  the  ra¬ 
tio  Eu/su  equal  to  1,200.  The  undrained  shear  modulus,  G  ,  is  one  third 
of  Eu  . 

55.  Data  for  the  computed  f-z  curves  are  given  in  Table  6.  The  as¬ 
sumed  pile  movements  in  the  first  column  were  used  for  the  first  three  methods 
shown.  To  compute  values  of  side  shear  for  the  method  by  Coyle  and  Reese 
(1966)  values  of  the  ratio  of  the  load  transfer  to  the  shear  strength,  ob¬ 
tained  by  using  "Curve  C"  in  Figure  10,  were  multiplied  by  648.  Equation  23 
was  used  with  f^ y  2  648  and  z  =  0.25  in.  to  compute  values  of  side  shear 
for  the  Vijayvergiya  (1977b)  method.  For  the  Kraft,  Ray,  and  Kawaga  (1981) 
method,  values  of  pile  movement  were  computed  using  Equation  27  for  selected 
values  of  side  shear  ranging  from  0.10  times  fmax  to  0.95  times  fmax  .  The 
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zone  of  influence,  rm  ,  was  computed  using  v  =  0.5  and  p  =  0 . 6 1 6  .  The 


parameter  4*  was  computed  assuming  R^  =  1  .  For  the  method  by  Heydinger  and 
O'Neill  (in  preparation),  values  of  K  =  7.395  and  m  =  5.672  were  deter¬ 
mined  using  Equations  30  and  31,  respectively,  with  L/D  =  75  .  An  average 
value  of  Eu  =  1,212,000  psf  was  computed  for  the  soil  along  the  total  length 
of  embedment.  Computed  values  of  side  shear  were  determined  for  the  assumed 
pile  movements  using  =  Eu/k  =  117,086  ,  D  =  12  in.  and  fmax  =  648  . 


Discussion  of  Predictive  Methods 


56.  The  method  by  which  the  four  predictive  methods  are  compared  by 
evaluating  them  according  to  their  ability  to  account  for  different  factors 
that  affect  the  shapes  of  f-z  curves.  Factors  affecting  f-z  curves  are 
concerned  with  pile  and  soil  properties.  The  pile  properties  that  have  been 
considered  are  length,  diameter,  L/D  ratio,  and  stiffness.  For  the  purpose 
of  the  discussion  on  pile  properties,  it  has  been  assumed  that  the  magnitude 
of  fmay  would  not  be  affected  by  pile  properties  as  reported  from  theoreti¬ 
cal  analyses  (Heydinger  and  O'Neill,  in  preparation).  Soil  properties  include 
undrained  shear  strength,  stiffness,  depth,  sensitivity,  and  stress  history. 


Pile  length 


57.  The  effects  of  pile  length  have  been  investigated.  Elasticity 
(Randolph  and  Wroth  1978)  or  finite  element  (Heydinger  and  O'Neill,  in  prep¬ 
aration)  solutions  indicate  that  soil  shears  caused  by  pile  loading  encompass 
larger  areas  radially  for  longer  piles.  The  integral  equation  of  the  dis¬ 
placement  in  the  shear  zones,  Equation  25,  indicates  that  displacements  for  a 
given  level  of  side  shear  would  be  greater  for  longer  piles.  This  implies 
that  f-z  curves  would  be  less  linear  and  that  displacements  required  to  mo¬ 
bilize  maximum  side  shear  would  be  greater  for  longer  piles. 

58.  Predicted  f-z  curves  would  be  affected  with  two  of  the  four  meth¬ 
ods.  According  to  the  Kraft,  Ray,  and  Kawaga  (1981)  method  rm  would  be 
greater  for  longer  piles.  Therefore,  predicted  values  of  pile  movement  using 
Equation  27  would  be  greater  for  longer  piles.  With  side  shear  plotted  as  the 
ordinate  and  pile  movement  as  the  abscissa,  f-z  curves  would  be  shifted  to¬ 
wards  the  right  as  pile  length  increased.  With  the  method  by  Heydinger  and 
O'Neill  (in  preparation),  the  parameter  K  increases  as  the  length  or  L/D 
increases  as  in  Equation  30;  therefore,  the  initial  slope  of  the  f-z  curves, 
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Efz  =  Eu/K  *  decreases.  The  parameter  m  ,  Equation  31,  decreases  as  L  in¬ 
creases  also,  and  the  f-z  curves  shift  to  the  right  as  m  decreases.  Both 
of  these  parameters  would  then  result  in  a  shift  to  right  of  the  f-z  curve 
as  pile  length  increases. 


Pile  diameter 


59.  Pile  diameter  affects  the  shape  of  f-z  curves  and  the  movement, 
z„  ,  required  to  mobilize  f__v  .  The  radial  zone  around  piles  where  soil  is 
remolded  during  installation  and  sheared  during  pile  loading  increases  as  pile 
diameter  increases.  Consequently,  larger  vertical  pile  displacements  are  re¬ 
quired  to  mobilize  side  shear  for  larger  diameter  piles.  These  findings  have 
been  verified  by  elasticity  solutions  and  finite  element  analyses.  Figure  14 
shows  the  results  of  finite  element  analyses  (Heydinger  and  O'Neill,  in  prep¬ 
aration)  where  diameters  were  varied.  By  increasing  diameters  (decreasing 
L/D  ratios),  the  f  versus  z/D  curves  shifted  to  the  left.  If  plots  of  f 


Figure  14.  Side  resistance  versus  z/D  ratio  from 
Finite  Element  Method  parametric  study,  SFBM  test 
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versus  z  were  made  for  the  curves,  the  initial  slopes  of  the  curves  would  be 
greater  for  the  larger  values  of  L/D  .  Thus,  there  is  no  unique  curve  for  f 
versus  z/D  or  for  f  versus  z  . 

60.  Three  of  the  methods  can  account  for  the  effect  of  diameter.  The 
method  by  Vijayvergiya  (1977b),  Equation  23,  would  result  in  a  shift  of  f-z 
curves  to  the  right  if  it  were  assumed  that  zc  increased  with  diameter.  For 
the  Kraft,  Ray,  and  Kawaga  (1981)  method,  with  Equation  24,  the  term  with 
rn/ro  ten8s  to  cause  zs  to  decrease  as  the  diameter  increases,  but  the 
overall  effect  is  for  zg  to  increase  as  diameter  increases  since  zg  is  di¬ 
rectly  proportional  to  rQ  .  According  to  the  method  by  Heydinger  and  O'Neill 
(in  preparation),  K  ,  Equation  30,  decreases  as  L/D  decreases  and  Efz 

=  Ey/K  increases.  The  shape  parameter,  Equation  30,  increases  as  L/D  de¬ 
creases.  Thus,  Efz  increases  and  m  increases  as  diameter  increases  (L/D 
decreases),  but  the  computed  side  shear  decreases  since  E^  is  multiplied 
by  z/D  in  Equation  29.  The  effect  of  increasing  the  diameter  then  is  to 
shift  f-z  curves  to  the  right.  Figure  15  indicates  the  variation  of  the 
computed  f-z  curves . 

L/D  ratio 

61.  An  additional  verification  to  the  effects  of  pile  length  and  diame¬ 
ter  on  the  initial  slopes  of  f-z  curves  is  obtained  in  terms  of  the  L/D 
ratio  from  the  work  of  others.  Theoretical  results  by  Baguelin  and  Frank 
(1980)  were  reported  giving  the  slope  of  the  initial  tangent  to  the  f  versus 
z/rQ  curves  as  a  function  of  the  ratio  G^k  ,  where  Gj^  is  the  initial 
shear  modulus.  The  parameter  k  is  shown  in  Figure  16  as  a  function  of 

L/D  .  The  figure  illustrates  that  the  slope  increases  as  L/D  increases. 
Equivalent  results  were  also  plotted  for  K  using  Equation  28  and  assuming 
that  E  =  3C-  .  Figure  16  indicates  that  there  is  close  agreement  between  re¬ 
sults  obtained  from  the  theoretical  studies  and  those  that  are  obtained  from 
the  method  by  Heydinger  and  O'Neill.  Similar  results  could  be  obtained  using 
the  method  by  Kraft,  Ray,  and  Kawaga  (1981). 

Pile  stiffness 

62.  Pile  stiffness,  usually  expressed  as  the  product  of  the  cross- 
sectional  area  and  Young  s  modulus  of  the  pile,  has  been  investigated 
(Heydinger  and  O'Neill,  in  preparation).  Values  of  the  stiffness  were  varied 
by  factors  of  10  and  1/10  for  the  finite  element  analyses  of  three  well- 
instrumented  pile-load  tests.  By  comparison  of  the  predicted  f-z  curves, 
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Figure  15.  Effect  of  pile  length  and  diameter  on  f-z  curves 

it  was  determined  that  pile  stiffness  does  not  nfluence  the  shapes  of  f-z 
curves  significantly.  However,  the  load-deformation  behavior  of  the  pile  top 
is  affected  by  its  stiffness  since  the  compressibility  changes.  The 
predictive  methods  consequently  do  not  include  pile  stiffness  as  a  variable. 
Undralned  shear  strength 

63.  For  the  purpose  of  this  comparison,  the  undrained  shear  strength  or 

fma„  is  considered  since  f_av  is  dependent  on  this  strength.  The  Coyle  and 
max  max 

Reese  (1966)  procedure  utilizes  curves  of  mobilized  side  resistance  dividea  by 
the  effective  undrained  shear  strength  which  is  equivalent  to  fmax  •  The 
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Figure  16.  Initial  slope  of  f-z  curves 
(Heydinger  and  O'Neill,  in  preparation) 

undrained  shear  strength  was  determined  from  unconfined  compression  tests.  The 
basis  for  the  three  recommended  curves  is  the  trends  that  were  observed  from 
field  and  laboratory  tests.  Vijayvergiya  (1977b)  proposed  a  parabolic  equa¬ 
tion  for  f/f_QV  as  a  function  of  z/z„  using  the  curves  recommended  by 
Coyle  and  Reese  (1966).  Hdwever,  efforts  to  normalize  f-z  curves  using 
f/fmay  did  not  result  in  strong  trends  for  the  initial  slopes  or  the  shapes 
of  the  f-z  curves  as  a  function  of  depth  or  soil  stiffness  for  the  three 
pile-load  tests  reported  by  Heydinger  and  O'Neill  (in  preparation).  In  any 
case,  computed  values  of  mobilized  side  resistance  are  directly  proportional 
to  fmav  for  the  methods  by  Coyle  and  Reese  (1966)  and  Vijayvergiya  (1977b). 

64.  The  other  two  methods  rely  on  the  stress  level  to  determine  the  de¬ 
gree  of  nonlinearity  of  the  f-z  curves.  The  stress  level  is  expressed  as 
the  ratio  f/fmax  •  For  the  Kraft,  Ray,  and  Kawaga  (1981)  procedure,  the  pa¬ 
rameter  which  is  proportional  to  the  stress  level,  occurs  in  Equation  27. 

For  the  method  by  Heydinger  and  O'Neill  (in  preparation),  the  term  Efz 


(z/D)/fnjax  in  the  denominator  of  Equation  29  is  equivalent  to  the  stress 
level.  For  both  methods  the  stress  level  influences  the  curvature  of  f-z 
curves.  Computed  values  of  z  or  f  are  not  directly  proportional  to  esti¬ 
mated  values  of  undrained  shear  strength  or  f,^  ,  so  the  procedures  are  not 
as  sensitive  to  errors  in  estimating  the  quantities. 

Soil  stiffness 

65.  Soil  stiffness  can  be  expressed  in  terms  of  Young's  modulus  or 
shear  modulus.  In  this  case,  undrained  moduli  were  used  to  represent  un¬ 
drained  pile  loading.  Methods  that  include  a  soil  modulus  correlate  the  in¬ 
teraction  behavior  at  the  pile-soil  interface  to  the  soil  stiffness.  There  is 
a  rational  basis  for  such  procedures  since  the  side  shear-deformation  behavior 
of  soil  near  the  pile  surface  is  related  to  the  stress-strain  properties  of 
soil. 

66.  Two  methods  rely  on  soil  moduli  to  compute  f-z  curves.  The  Kraft, 
Ray,  and  Kawaga  (1981)  solution  incorporates  a  hyperbolic  expression  to  repre¬ 
sent  the  average  shear  modulus,  which  is  dependent  on  the  level  of  stress. 

Since  z  is  inversely  proportional  to  the  initial  shear  modulus  in  Equa¬ 
tion  27,  error  involved  in  estimating  GQ  is  directly  evident  throughout  the 
computations.  The  general  equation,  Equation  29,  has  both  theoretical  and  em¬ 
pirical  basis  for  computing  f-z  curves  as  a  function  of  Eu  .  With  this 
method,  errors  in  estimating  Ey  result  in  errors  for  Efz  and  the  shape  pa¬ 
rameter  m  .  However,  the  error  diminishes  as  the  stress  level  increases. 

Figure  17  illustrates  the  effects  of  varying  soil  stiffness  for  the  two  methods. 
Soil  depth 

67.  Soil  depth  affects  soil  confining  pressures.  Coyle  and  Reese 
(1966)  showed  that  load  transfer  increased  with  depth.  They  attributed  the 
effect  of  depth  to  the  fact  that  there  are  vibrations  during  driving  which 
open  small  spaces  between  the  pile  and  the  soil  which  may  not  close  near  the 
ground  surface  after  driving.  Therefore,  three  curves  were  recommended  de¬ 
pending  on  depth. 

68.  According  to  theoretical  approaches,  side  shear  is  dependent  on 

soil  shear  strength  and  confinement,  and  side  shear-deformation  behavior  is 
dependent  on  the  soil  stiffness.  When  using  Equation  27,  the  computed  pile 
movements  decrease  when  or  the  soil  modulus  increase.  Similarly,  the 

computed  mobilized  side  shear  increases  when  fmax  or  the  soil  modulus  in 
Equation  29  is  increased.  Thus,  assuming  one  or  both  of  the  factors  increase 


with  depth,  the  f-z  curves  would  shift  upward  as  could  be  expected.  The 
methods  cannot  account  for  spaces  between  pile  and  soil  if  they  exist.  Side 
shear  should  be  ignored  where  the  spaces  exist. 

69.  The  measured  results  from  three  pile-load  tests  were  used  to  de¬ 
velop  the  parameters  for  the  general  equation.  Typically  the  f-z  curves 
were  shifted  upward  with  increases  in  depth.  However,  when  attempting  to  cor¬ 
relate  K  and  m  with  depth,  it  was  determined  that  no  definite  correlation 
between  the  parameters  could  be  made  with  depth.  Therefore,  the  same  values 
of  K  and  M  should  be  used  for  a  given  pile  unless,  for  the  case  of  m  , 
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there  are  strata  with  significant  differences  in  the  undrained  soil  modulus. 
Soil  sensitivity  and  remolding 


70.  Soil  sensitivity  affects  the  properties  of  soils  when  they  are 
remolded  during  pile  installation.  It  has  been  reported  (Coyle  and  Reese 
1966)  that  soft,  insensitive  clays  can  have  higher  undrained  shear  strengths 
after  the  soil  reconsolidates  around  piles.  On  the  other  hand,  triaxial  test¬ 
ing  on  sensitive  clays  (Roy  et  al.  1981)  have  shown  that  the  strengths  and 
moduli  of  remolded  clays  were  reduced  by  30  and  65  percent,  respectively,  due 
to  remolding.  However,  it  has  not  been  shown  conclusively  that  soil  prop¬ 
erties  after  reconsolidation  would  be  significantly  reduced. 

Stress  history 

71.  Stress  history  is  represented  by  the  overconsolidation  ratio. 
Overconsolidation  can  occur  as  the  result  of  glaciation,  sediment  erosion, 
dessication,  and  because  of  secondary  consolidation  effects.  Profiles  of  nor¬ 
mally  consolidated  soils  consist  of  soils  with  relatively  homogeneous  soil 
properties  with  depth  or  gradually  changing  with  depth.  Profiles  of  these 
soils  generally  exhibit  a  marked  decrease  in  overconsolidation  effects  with 
depth.  Methods  by  Coyle  and  Reese  (1966)  and  Vijayvergiya  (1977b)  relying 
principally  on  s,,  or  fm_„  to  compute  f-z  curves  do  not  account  for  over- 
consolidation  effects  adequately.  The  other  two  methods  make  use  of  soil  mod¬ 
uli  which  account  for  overconsolidation  effects.  Figures  18  and  19  illustrate 
predicted  f-z  curves  by  the  four  methods  for  normally  consolidated  and  over¬ 
consolidated  clays,  respectively.  The  data  listed  in  Table  6  were  used  to 
plot  the  curves  in  Figure  18. 

Comparison  of  Predictive  Methods 

72.  The  four  instrumented  pile-load  tests  that  were  used  to  evaluate 
the  proposed  methods  for  predicting  side  shear  were  also  used  for  comparing 
f-z  curves.  Representative  f-z  curves  were  selected  for  each  of  the  four 
tests.  Soil  properties  were  estimated  from  the  reports  of  the  pile-load 
tests.  For  the  first  three  tests,  the  soil  properties  were  those  interpreted 
for  the  finite  element  analyses  that  were  used  to  develop  the  general  equation 
(Heydinger  and  O'Neill,  in  preparation).  The  a-method  recommended  by  Dennis 
and  Olson  (1983)  was  used  to  compute  fmax  for  the  comparisons. 


UHCC  test  (O’Neill, 

Hawkins,  and  Mahar  1981) 

73.  Computed  f-z  curves  were  obtained  for  two  depths,  21  and  36  ft. 
Values  of  Eu  corresponding  to  values  of  Eu/su  equal  to  1,000  and  900,  re¬ 
spectively,  for  the  two  depths.  The  values  of  Eu  are  higher  than  values  in¬ 
terpreted  from  triaxial  or  pressuremeter  tests  and  are  lower  than  values  in¬ 
terpreted  from  crosshole  shear  tests.  The  computed  results  were  compared  to 
measured  results  at  nearby  depths.  The  comparisons  are  shown  in  Figure  20. 
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Figure  19.  Predicted  f-z  curves  for  overconsolidated  clay 


SFBM  test  (Kirby  and 
Roussel  1979  and  1980) 


74.  The  computed  and  measured  results  are  for  depths  of  20  and  35  ft. 
Values  of  Eu/su  equal  to  600  were  used  to  estimate  Eu  ,  as  recommended  by 
Kirby  and  Roussel  (1979  and  1980).  They  are  representative  of  moduli  obtained 
from  triaxial  tests  at  low  strains.  The  comparisons  are  given  in  Figure  21. 
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St.  Alban  test  (Blanchet, 
Tavenas,  and  Garneau  1980; 


Konrad  1977;  Rov  et  al.  1981) 

75.  Computed  f-z  curves  were  compared  to  a  curve  computed  from  the 
average  measured  side  friction  since  strain-gage  measurements  were  not  made. 

A  ratio  of  equal  to  450  was  used  as  recommended  (Roy  et  al.  1981)  for 

remolded  soil  on  triaxial  testing.  Values  of  Eu/su  equal  to  900  were  ob¬ 
tained  from  undisturbed  or  aged  samples.  The  curves  are  given  in  Figure  22. 

EABPL  test  (US  Army  Engineer 
District,  Mew  Orleans  1977; 

Shilstone  Testing  Laboratory  1976) 

76.  Computed  and  measured  f-z  curves  were  obtained  for  depths  of  17 
and  49  ft.  Eu  was  estimated  using  Eu/su  equal  to  600.  A  wide  range  of 
values  of  the  ratio  Eu/su  was  obtained  from  triaxial  testing,  thereby  re¬ 
quiring  a  reasonable  estimate.  The  comparisons  are  shown  in  Figure  23. 

Soil-Structure  Interaction  Program 

77.  The  purpose  of  this  section  is  to  illustrate  the  effects  that  the 
shapes  of  f-z  curves  have  on  the  predicted  load-displacement  behavior  of  the 
pile  top.  Of  particular  interest,  were  the  effects  of  using  nonlinear  f-z 
curves  of  varying  degrees  of  nonlinearity.  This  aspect  was  investigated  for 
soils  with  properties  representative  of  soft  and  stiff  clays.  The  effects  of 
varying  the  initial  slopes  of  f-z  curves  were  also  investigated  to  illus¬ 
trate  the  sensitivity  of  pile-top  behavior  on  the  initial  slope.  The  results 
indicate  the  importance  of  accurate  determination  of  f-z  curves. 

78.  For  this  investigation,  the  axially  loaded  pile  analyses  were  con¬ 
ducted  using  the  program  titled  PX4C3  (Radhakrishnan  and  Parker  1975),  which 
was  developed  at  the  University  of  Texas  at  Austin.  The  program  requires 
input  for  the  pile  tip  load-deformation  behavior,  referred  to  as  the  Q-z 
curve,  and  for  a  number  of  f-z  curves.  In  order  to  investigate  the  effects 
of  the  f-z  curve  input,  Q-z  curve  data  were  input  and  then  the  f-z 
curves  were  varied.  The  criteria  used  to  determine  the  Q-z  data  were  given 
in  the  documentation  for  the  program  (Radhakrishnan  and  Parker  1975)  for  the 
program  for  soils  with  no  available  stress-strain  curve.  A  50-ft-long  pile 
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with  an  outside  diameter  of  1  ft  and  stiffness  AE  =  4.37  x  10  was  used 
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Figure  22.  Measured  and  computed  f-z  curves 
for  St.  Alban  test 


throughout,  f-z  curves  were  computed  for  five  depths  at  10-ft  intervals  be 
ginning  at  a  depth  of  5  ft. 

79.  The  variation  of  the  f-z  curves  was  achieved  using  the  criteria 
by  Heydinger  and  O'Neill  (in  preparation).  With  this  method,  the  single 
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parameter  m  can  be  used  conveniently  to  vary  the  shapes  of  the  f-z  curves, 
and  the  parameter  Efz  can  be  used  to  vary  the  initial  slope.  Similar  re¬ 
sults  can  be  obtained  by  varying  the  parameters  for  the  method  by  Kraft,  Ray, 
and  Kawaga  (1981). 

Case  1 — variation  of 

f-z  curve  shapes,  soft  soils 

80.  The  behavior  of  a  pile  in  a  soft  soil  was  modeled  using  represen¬ 
tative  soil  parameters  for  a  saturated  clay.  A  saturated  unit  weight  of 
102.4  pcf  was  used  to  compute  the  effective  vertical  overburden  pressure  dis¬ 
tribution.  The  undrained  shear  strength  was  computed  at  various  depths  using 
the  ratio  of  the  undrained  shear  strength  to  the  effective  vertical  overburden 
pressure  equal  to  0.35.  The  undrained  soil  modulus  was  computed  using  the  ra¬ 
tio  Eu/su  equal  to  800  and  two  types  of  soil  profiles  were  investigated  for 
this  case.  Profile  1  consists  of  a  soil  formation  with  soil  strengths  and 
moduli  increasing  with  depth  as  determined  by  the  input  parameters  described 
above.  Profile  2  consists  of  a  single,  homogeneous  layer  of  soil  having  the 
parameters  determined  for  a  depth  of  25  ft  for  Profile  1 . 

81.  The  input  data  for  the  Q-z  curve  for  Profile  1  was  obtained  by 
using  the  computed  soil  modulus  for  a  depth  of  50  ft.  For  this  case  then, 

Eu  =  560,000  psf  was  computed.  The  corresponding  stress-strain  curve  and  the 
Q-z  data  were  computed  using  the  recommended  criteria  ( Radhakrishnan  and 
Parker  1975).  The  same  Q-z  data  were  used  for  Profile  2. 

82.  Side  shear-displacement  curves  were  input  for  the  five  depths  indi¬ 
cated  previously  for  Profiles  1  and  2.  According  to  Equation  30,  the  parame¬ 
ter  K  is  6.34  for  L/D  =  50  for  all  depths  of  Profiles  1  and  2.  The  ini¬ 
tial  slope  of  the  line  that  is  tangent  to  the  curve  of  f  versus  z  was  com¬ 
puted  using  Eu/k  .  Thus,  for  example,  the  average  undrained  soil  modulus  at 
the  25-ft  depth,  Eu  =  280,000  ,  was  used  to  compute  the  initial  slope  to 

the  f-z  curves,  E^z  -  44,164  .  The  shape  parameter,  m  =  3.14  ,  was  com¬ 
puted  using  Equation  31  and  the  average  value  of  Eu  .  The  f-z  curves  were 
computed  using  fmax  =  su  in  Equation  29.  Figure  24  shows  the  computed 
f-z  curves  for  different  values  for  m  for  a  depth  of  25  ft.  It  indicates 
the  variation  of  the  curves  with  m  and  the  curve  with  m  =  10  approaches  a 
bilinear  curve.  Similar  sets  of  curves  were  obtained  for  the  other  depths. 

83.  The  results  of  the  axially  loaded  pile  analyses  are  shown  in  Fig¬ 
ure  25.  For  the  analyses,  sets  of  f-z  curves  were  input  for  each  of  the 
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Figure  24.  Variation  of  f-z  curves  with  m 


values  of  m  shown  in  Figure  24.  The  pile  head  deformation  variation  for 
Profile  1  was  0.15,  0.08,  0.07  in.  for  m  =  1,0  ,  3.14,  and  10,  respectively, 
for  the  working  loads  which  were  taken  as  one  half  the  ultimate  loads.  Also, 
shown  in  Figure  25  is  the  variation  of  the  pile  head  deformations  for  Pro¬ 
file  2.  The  deformations  at  the  working  loads  are  0.10,  0.07,  and  0.06. 

Case  2 — variation  of 

f-z  curve  shapes,  stiff  soil 

84.  A  homogeneous  soil  profile  was  used  for  the  stiff  soil  analyses 
similar  to  the  method  that  was  used  for  Profile  2  for  the  soft  soil  analyses. 
A  saturated  unit  weight  of  122.4  pcf  was  used  to  compute  the  vertical  effec¬ 
tive  overburden  pressure  for  a  depth  of  25  ft,  equal  to  1,500  psf.  The  ratio 
of  the  undrained  shear  strength  to  the  effective  vertical  pressure  equal  to 
0.8  was  used  to  compute  su  =  1,200  psf.  An  undrained  soil  modulus,  Eu 

=  1,440,000  psf,  was  computed  using  Eu/su  =  1,200  . 

85.  The  data  for  the  Q-z  and  f-z  curves  were  computed  as  described 
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Figure  25.  Pile  top  load  versus  pile  top 
movement,  Case  1,  Profiles  1  and  2 


previously.  The  stress-strain  curve  was  computed  by  using  the  undrained  soil 
modulus.  The  parameter  K  =  6.34  remains  unchanged.  The  initial  slope  is 
Efz  =  Eu  /K  =  227,129  .  The  computed  shape  parameter  is  7.61.  The  side  shear- 
displacement  curves  were  computed  using  fmax  =  0-8,  su  =  960  . 

86.  Figure  26  shows  the  results  of  the  pile  analyses  for  f-z  curves 
with  m  =  7.61  and  for  bilinear  functions.  The  pile  head  deformation  varied 
from  0.07  in.  for  the  m  =  7.61  to  0.06  in.  for  the  bilinear  function  for  the 
working  load. 

Case  3 — variation  of  ini- 
tial  slope  of  f-z  curve 

87.  The  analyses  with  the  variation  of  the  initial  slopes  of  the  f-z 
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Figure  26.  Pile  top  load  versus  pile 
top  movement,  Cases  1  and  2 

curves,  Efz  ,  were  conducted  for  the  previous  analyses  of  Case  2  by  using  m 
=  7.61  .  The  initial  slopes  were  varied  by  factors  of  0.5  and  2  and  the  re¬ 
sulting  f-z  curves  were  input.  The  results  of  the  analyses  are  shown  in 
Figure  26.  The  pile-top  deformation  for  2Efz  is  0.05  in.  and  the  defor¬ 
mation  for  (1/2)Efz  is  0.09  in.  The  differences  in  the  pile  deformations  at 
the  design  load  would  have  been  larger  had  bilinear  curves  been  used.  For 
these  three  cases,  the  magnitude  of  the  variations  of  the  pile-top  deforma¬ 
tions  are  not  large,  but  they  are  significant  if  the  percent  differences  are 
considered.  The  differences  would  be  larger  for  larger  diameter  piles  or 
longer  piles,  or  if  loads  are  to  exceed  the  working  loads. 


Condition  and  Recommendations 


88.  The  methods  by  Coyle  and  Reese  (1966)  and  Vijayvergiya  (1977b)  have 
many  limitations.  They  do  not  account  for  factors  which  affect  the  shapes  of 
f-z  curves,  such  as  soil  stiffness  and  L/D  ratio.  Additionally  computed 
values  of  side  friction  are  directly  proportional  to  fMX  ,  so  the  shapes  of 
the  curves  are  affected  by  errors  in  computing  f  .  Therefore,  the  two 
methods  are  not  recommended  for  use. 

89-  The  methods  by  Kraft,  Ray,  and  Kawaga  (1981)  and  Heydinger  and 
O'Neill  (in  preparation)  are  comparable.  They  both  have  theoretical  basis  and 
have  been  verified  with  comparisons  from  field  tests.  The  expression  for  the 
Heydinger  and  O'Neill  method  comes  from  a  general  equation.  Kraft,  Ray,  and 
Kawaga  (1981)  use  a  hyperbolic  equation  to  represent  the  change  of  the  shear 
modulus,  which  results  in  f-z  curves  with  shapes  that  are  hyperbolic.  Both 
methods  account  for  factors  which  affect  the  shapes  of  f-z  curves.  They  are 
both  sensitive  to  the  stress-strain  modulus  that  is  input. 

90.  The  recommendation  is  to  use  either  the  method  by  Kraft,  Ray,  and 
Kawaga  (1981)  or  the  method  by  Heydinger  and  O'Neill  (in  preparation).  The 
computed  curves  are  very  similar.  It  has  been  shown  (Heydinger  and  O'Neill, 
in  preparation;  and  Kraft,  Ray,  and  Kawaga  1981)  that  the  two  methods  can  be 
used  to  predict  f-z  curves  very  accurately  if  representative  values  of  Eu 
are  known.  Further  verifications  should  be  made  using  data  from  other  tests 
wherever  possible  in  order  to  come  to  a  better  understanding  of  the  methods 


PART  IV:  ENGINEERING  PROPERTIES  OF  SOILS 


Introduction 


91.  The  purpose  of  this  section  is  to  describe  research  on  engineering 
properties  of  soil.  Soil  properties  can  be  estimated  by  using  correlations 
with  soil  plasticity  or  normalized  soil  testing.  The  concept  of  normalized 
soil  properties  has  been  found  to  be  valid  for  many  soils.  The  following 
findings  were  obtained  from  research  on  soil  properties  for  general  geotechni¬ 
cal  purposes  and  not  specifically  for  use  in  design  of  piles. 

Correlations  With  Plasticity 

Correlations  for  ♦' 

92.  Correlations  for  the  effective  angle  of  internal  friction  have 
been  proposed  by  a  number  of  investigators  for  normally  consolidated  clays. 
Figure  27,  obtained  from  NAVFAC  DM-7  (1979),  shows  the  results  reported  by 
BJerrum  and  Simons  (I960).  The  true  angle  of  internal  friction,  $r  ,  is 


determined  from  the  deviator  stress  from  samples  with  the  same  water  contents. 
The  curve  for  the  effective  friction  angle,  ♦'  ,  was  obtained  from  drained 
tests  and  from  undrained  tests  with  pore  pressure  measurements.  For  the  un¬ 
drained  tests,  was  determined  using  the  maximum  value  of  the  ratio  of  the 

maximum  stress  to  the  minimum  stress,  /o^  >  as  the  failure  criterion. 

93.  Mayne  (1980)  reported  the  following  equation  for  for  both  nor¬ 

mally  and  overconsolidated  clays  based  on  the  plasticity  index,  PI  ,  and  the 
liquid  limit,  LL 

sin  ♦'  =  0.656  -  0.409  jj  (32) 

The  correlation  coefficient  for  Equation  32  is  0.583.  Figure  28  shows  the  to¬ 
tal  friction  angle  obtained  from  consolidated-undrained  triaxial  tests,  $  , 

for  marine  soils.  Based  on  these  findings,  it  is  apparent  that  the  effective 
friction  angle  correlates  poorly  with  index  properties. 


Plasticity  index,  PI,  l 

Figure  28.  Correlation  for  ^cu( Vijayvergiya  1977b) 
Correlations  for  KQ 


94.  The  ratio  of  the  horizontal  effective  stress  to  the  vertical  ef¬ 
fective  stress,  Kq  ,  is  dependent  on  soil  plasticity  and  stress  history. 
Figure  29,  from  Ladd  et  al.  (1977),  indicates  a  general  trend  for  KQ  as  a 
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function  of  plasticity  index  and  ♦  '  for  normally  consolidated  soils.  The 
following  equation  by  Alpan  (1967)  gives  KQ  in  terms  of  plasticity  for  nor¬ 
mally  consolidated  clays. 

K  =  0.19  ♦  0.233  log  PI(%)  (33) 

o 

95.  An  expression  for  the  ratio  of  KQ  for  overconsolidated  soils, 

(K_ )  ,  to  the  value  for  the  soil  if  it  was  normally  consolidated,  ( K  )  , 

°  oc  nc 

is  given  in  terms  of  the  OCR  . 

=  0CRn  (3*0 

Figure  30  contains  a  curve  for  the  parameter  n  as  a  function  of  PI  for 
undisturbed  soils  (Vesic  1972).  Figure  31  was  proposed  by  Brooker  and  Ireland 
(1965).  Vijayvergiya  (1977a)  computed  Figure  32  using  Figure  27  to  compute 
+  and  the  following  equation. 

Kq  =  (1  -  sin  ^KOCR)172  (35) 


Figure  30.  Exponent  n  versus 
PI  (Vesic  1977) 


Plasticity  index,  PI,  7. 

Figure  31.  KQ  as  a  function  of  PI  and  OCR 
(Brooker  and  Ireland  1965) 

The  figures  and  equations  indicate  that  KQ  increases  with  increasing  plas¬ 
ticity  for  both  normally  and  overconsolidated  soils,  and  that  KQ  increases 
with  Increasing  OCR  . 

96.  In  addition  to  the  previous  recommendations,  a  statistical  study  by 
Mayne  and  Kulhawy  (1982)  has  also  been  reported.  The  best  relationship  for 
KQ  for  normally  consolidated  clays  was  given  as 


987  sin 


The  correlation  coefficient  for  Equation  36  is  0.854.  The  authors  also  at¬ 
tempted  to  correlate  (K  )  with  plasticity  and  other  index  properties  but 

0  nc 

concluded  that  no  good  relationships  existed.  Mayne  and  Kulhawy  (1982)  also 
recommended  KQ  for  overconsolidated  clays. 

60  =  (1  -  sin  ♦')0CRsin  ♦  ’  (31 
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Figure  32. 


KQ  as  a  function  of  PI  and  OCR 


(Vijayvergiya  1977a) 


Correlations  for  un¬ 
drained  shear  strength 


97.  Relationships  have  been  proposed  for  the  undrained  shear  strength 
normalized  by  the  vertical  effective  stress,  su/°v  »  as  a  function  of  plas¬ 
ticity.  Using  the  results  of  vane  shear  and  unconfined  compression  tests, 
BJerrum  and  Simons  (I960)  have  shown  in  Figure  33  that  the  ratio  su/p  , 

(where  p  =  a^),  increases  as  the  PI  increases  for  normally  consolidated  ma¬ 
rine  clays.  They  also  showed  in  Figure  34  that  the  ratio  decreases  as  the  li¬ 
quidity  index  increases,  and  stated  that  values  of  the  ratio  would  be  higher 
for  freshwater  soils.  Figure  35,  from  NaKase  and  Kamei  (1983)  indicates  re¬ 
lationships  for  the  ratios  for  compression  and  extension  tests  for  freshwater 
soils.  Curves  were  determined  for  normally  consolidated  freshwater  soils  from 
the  results  of  triaxial  compression  (TC)  and  direct  simple  shear  tests  (DSS) 
as  shown  in  Figure  36.  Mayne  (1980)  proposed  the  following  relationship  for 
normally  consolidated  clays. 


le*  V  «■  ■» 

l,-»  V* 

I  fS  i.'  *  *  *. 

iv«\v,  iv 


I  \\v 

»•  » 


1  iV\vav«*. 

I&V-- 


i.i 


te 

\$m 

RvM'w 


lil 


wsww 
vW »> 


</J  |o.  0.20 


■Calculated  from  trlaxlal  tests 
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Figure  33.  su/p  versus  PI 


(BJerrum  and  Simons  I960) 


Liquidity  index. 


w.  -  w 
L  p 


Figure  34.  su/p  versus  LI 
(BJerrum  and  Simons  I960) 


w 

-T  =  0.642  sin  <D'  +  0.031 


The  correlation  coefficient  for  Equation  38  is  0.722.  The  trend  is  for  su/°y 


to  increase  with  increasing  plasticity  and  to  decrease  as  the  water  content, 
i.e.,  LI  increases. 


Normalized  Soil  Testinc 


98.  The  results  presented  for  normalized  testing  are  based  on  the 
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Plasticity  index,  PI,  % 
Figure  35.  su/p  versus  PI 

(NaKase  and  Kamei  1983) 


normalized  soils  properties  (NSP)  concept  and  the  stress  history  and  normal¬ 
ized  soil  engineering  properties  (SHANSEP)  method  of  design  (Ladd  and  Foott 
1974).  The  hypothesis  is  that  the  effects  of  soil  overconsolidation  can  be 
duplicated  in  a  laboratory  by  consolidating  samples  to  pressures  greater  than 
any  previously  encountered.  After  the  samples  are  allowed  to  swell  back  to 
lower  pressures,  they  are  then  sheared  to  failure.  The  testing  procedure  is 
an  attempt  to  overcome  the  effects  of  soil  disturbance  due  to  sampling.  Soils 
that  are  sensitive  or  that  are  naturally  cemented  due  to  thixotropic  effects 
would  not  behave  according  to  normalized  behavior. 

Normalized  undrained  shear  strength 


99.  The  normalized  undrained  shear  strength  is  dependent  on  the  OCR 
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Figure  36.  su/p  versus  PI  for  normally 
consolidated  soil  (Gardner  1977) 

of  the  soil  (Simons  I960,  Skempton  1957).  By  consolidating  samples  and  then 
allowing  them  to  swell  back,  soils  with  known  OCR's  are  obtained.  The  OCR 
is  then  defined  as  the  ratio  of  maximum  vertical  consolidation  pressure, 

,  to  the  vertical  consolidation  pressure  after  swelling,  .  Typical 
curves  of  su/o^  versus  the  log  of  OCR  ,  shown  in  Figure  37,  indicate  that 
the  ratio  increases  similarly  for  soils  of  different  plasticity  (Vijayvergiya 
1977a).  There  is,  however,  no  definite  relationship  between  PI  and  the  ra- 

versus  the  log  of  OCR  ,  Vijayvergiya 


tio.  By  plotting  su/°y  su/ov 

oc  nc 


(1977a)  obtained  a  range  of  values  that  can  be  used  to  estimate  su/o^ 


for 


over consol Ida ted  soils  as  shown  in  Figure  38. 

100.  Expressions  for  su/o^  have  been  obtained  for  overconsolidated 
soils.  Gardner  (1977)  recommended  the  following  equation  for  CIU  triaxial 
testing  on  silty  clays  with  PI  =  16  ±  15  . 
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pm 


(39) 


■p  =  (0.111  ±  0.06)0CR 


Researchers  in  Japan  (Hanzawa  1977,  Hanzawa  and  Kish Ida  1982,  NaKase  and  Kamei 
1983)  recommended  the  following  equation  obtained  from  KQ  consolidation 
tests  (CK0U). 


s 

~r  =  (0.355)0CR  (110) 

v 


Overconsol idation  ratio,  OCR 

Figure  37.  3u/°y  versus  OCR 

(Vijayvergiya  1977a) 

very  low  strain, 
problems . 


101.  Researchers  have  reported 
that  soil  consolidated  to  high  pres¬ 
sures  in  the  laboratory  experiences  an 
increase  in  strength  with  time  (BJerrum 
and  Lo  1963,  Leroueil  et  al.  1979, 
Tavenas  and  Leroueil  1977).  The  expla¬ 
nation  is  that  soils  subjected  to  thix¬ 
otropic  or  secondary  consolidation  ef¬ 
fects  will  gain  strength  if  permitted 
to  age.  Reductions  of  shear  strength 
of  30  percent  for  young  clays,  as  com¬ 
pared  to  old  clays,  have  been  docu¬ 
mented.  The  implication  is  that  the 
normalized  testing  procedure  would  re¬ 
sult  in  reductions  in  shear  strength. 
Normalized  undrained  modulus 

102.  Undrained  soil  moduli  have 
been  normalized  by  the  undrained  shear 
strength  and  the  vertical  effective 
consolidating  pressures.  The  ratios 
are  dependent  on  the  OCR  of  the  soil 
and  on  plasticity.  The  ratios  also 
depend  on  the  level  of  shear  stress 
that  is  used  to  compute  the  undrained 

The  initial  modulus,  taken  at 
f-z  curves  and  other 


modulus, 
is  recommended  for  predicting 
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OCR  (Vijayvergiya  1977a) 


versus 


103.  The  ratio  Eu/su  has  been  found  to  be  very  sensitive  to  OCR  and 
plasticity.  The  ratio  decreases  appreciably  as  both  OCR  and  plasticity 
increase.  Values  ranging  from  1,200  to  1,500  have  been  measured  in  laboratory 
testing  for  soils  with  PI  less  than  30  percent  (D1 Appolonia,  Poulos,  and 
Ladd  1971).  The  following  equation,  based  on  analysis  of  a  number  of  tests, 
was  recommended  from  the  initial  shear  modulus,  GQ  (Hara  et  al.  1974) 

Gq  =  487su0'928  (kg/cm2)  (41) 

An  approximation  is  GQ/su  =  500  .  Values  of  Eu/su  as  low  as  100  have 
been  reported  for  highly  plastic  and  organic  soils  ( D * Appolonia,  Poulos,  and 

Ladd  1971). 
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104.  The  ratio  Eu/V  is  not  very  sensitive  to  overconsolidation  and 
plasticity.  Amarasinghe  and  Parry  (1975)  reported  that  the  ratio  did  not 
change  much.  Yudbir  and  Varadarajan  (1974)  reported  that  the  ratio  increases 
once  the  OCR  exceeds  a  value  of  100.  The  trend  for  ratios  of  Eu/o^  is 
that  they  decrease  somewhat  with  increasing  plasticity.  Typical  values  of  the 
ratio  are  in  the  range  of  100  to  400. 

Recomnendat ions 

105.  Some  specific  recommendations  for  estimating  the  engineering  prop¬ 
erties  of  soils  are  given,  followed  by  examples  of  their  use.  As  can  be  seen 
in  the  previous  discussion,  the  soil  properties  do  not  correlate  very  well. 
Therefore,  it  is  emphasized  that  the  recommendations  are  intended  only  to  be 
used  for  initial  estimates  or  as  an  aid  in  interpreting  soils  data.  They 
should  be  used  for  final  design  only  if  it  is  not  practical  to  obtain  quality 
triaxial  testing  results.  Thus,  it  is  still  necessary  to  conduct  extensive 
laboratory  tests  for  large  projects. 

106.  The  recommendations  are  based  on  the  reported  statistical  analyses 
and  on  the  results  of  normalized  testing.  The  recommendation  for  estimating 
the  effective  angle  of  internal  friction  is  to  use  Equation  32.  This  correla¬ 
tion,  although  not  very  good,  is  the  best  correlation  available.  The  two  cor¬ 
relations  for  KQ  that  can  be  used  for  normally  consolidated  soils  are  given 
in  Equations  33  and  36,  depending  on  the  information  that  is  available.  Equa¬ 
tion  35  or  37  can  be  used  to  estimate  KQ  for  overcj.isolidated  soils. 

107.  The  undrained  shear  strength  is  estimated  from  normalized  parame¬ 
ters.  For  normally  consolidated  soils,  the  ratio  su/°y  is  obtained  from 

Equation  38,  or  s  /o'  =  0.35  can  be  used.  The  ratio  s  /o'  can  be  obtained 
M  u  v  u  v 

for  overconsolidated  soils  using  an  equation  that  is  similar  to  Equation  39, 
where  the  quantity  0.41  +  0.06  is  equal  to  the  value  of  the  ratio  for  nor¬ 
mally  consolidated  soils.  By  plotting,  for  example,  the  curves  in  Figure  37 
using  the  logarithm  of  s^/o^  versus  the  logarithm  of  OCR  ,  a  set  of  lines 
that  are  approximately  parallel  is  obtained.  The  slope  of  the  lines  is  the  ex¬ 
ponent  of  OCR  in  Equation  39.  The  average  slope  of  the  lines  in  Figure  37 
is  approximately  0.73.  Therefore,  the  following  equation  is  recommended 


The  undrained  shear  strength  could  possibly  be  overestimated  using  Equation  40 

108.  The  undrained  soil  modulus  is  estimated  using  Eu/su  or  Eu/o^  * 

where  Eu  is  the  initial  modulus  or  the  slope  of  the  initial  tangent  to  the 
stress-strain  curves.  The  OCR  and  soil  plasticity  both  have  an  effect  on  the 
values  of  the  two  ratios,  with  the  ratios  tending  to  decrease  as  the  OCR  and 
the  PI  increase.  There  appears  to  be  no  consistent  trend;  therefore,  no  rela¬ 
tionships  are  given  in  this  report.  The  range  of  values  of  Eu/su  is  typi¬ 
cally  600  to  1,500.  The  range  of  values  of  Eu/oy  is  approximately  100  to 
400.  The  previous  analyses  with  f-z  curves  indicate  that  more  accurate 

f-z  curves  are  obtained  if  the  higher  values  in  the  ranges  of  values  are  used 

for  insensitive  soils,  and  if  lower  values  are  used  for  slightly  sensitive  and 
sensitive  soils. 

Example  of  predic- 
tion  of  soil  properties 

109.  The  first  example  of  prediction  of  soil  properties  was  made  for 

soil  at  the  SFBM  site  (Kirby  and  Rousell  1979  and  1980).  The  soil  is  a 
slightly  sensitive,  normally  consolidated  clay  of  high  plasticity.  The  com¬ 
puted  vertical  effective  stress  for  a  depth  of  20  ft  is  944  psf.  The  effec¬ 
tive  angle  of  internal  friction  was  computed  using  Equation  32  with  PI  =  37 

and  LL  =  93  .  The  computed  value  is  29.5  deg,  as  compared  to  an  average 

value  of  34  deg  determined  from  triaxial  testing.  Estimated  values  of  KQ 
equal  to  0.56  and  0.51  were  obtained  using  Equations  33  and  36,  respectively, 
using  29.5  deg  for  the  effective  friction  angle.  The  investigators  for  the 
pile  test  estimated  K0  =  0.49  .  The  ratio  of  su/cT  was  estimated  to  be 
0.346  using  Equation  38  as  compared  to  0.35  which  was  determined  from  triaxial 
testing.  The  ratio  Eu/*u  =  600  that  was  recommended  by  the  investigators  of 
the  site  corresponds  to  Eu/oy  =  208  .  Using  the  estimated  soil  properties, 
predicted  f-z  curves  similar  to  those  shown  in  Figure  21  would  be  obtained. 

110.  The  second  example  of  predictions  was  for  the  UHCC  site  (O'Neill, 
Hawkins,  and  Mahar  1981).  The  soil  at  the  site  is  an  insensitive,  highly 
overconsol idated  clay  of  high  plasticity.  The  effective  vertical  stress  at  a 
depth  of  20  ft  was  computed  to  be  1,912  psf.  The  effective  angle  of  internal 
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friction  equal  to  24.8  deg  was  computed  using  Equation  32  with  PI  =  40  and 
LL  =  68  .  An  effective  angle  equal  to  24  deg  was  determined  from  triaxial 
tests.  Kq  equal  to  1.42  and  1.23  was  computed  using  Equations  35  and  37. 
Values  of  K0  as  low  as  1  were  obtained  from  laboratory  tests,  and  values  of 
K0  as  high  as  2.5  were  obtained  from  pressuremeter  tests  at  the  site. 

111.  In  order  to  obtain  su/<jy  for  the  overconsolidated  soil, 
au/ov  =  wa3  estimated  for  a  normally  consolidated  soil  using  Equa¬ 

tion  38.  su/Oy  =  1-05  was  then  computed  using  Equation  42.  The  estimated 
undrained  shear  strength  would  then  be  2,008  psf  as  compared  to  1,700  which 
was  determined  from  normalized  testing  at  the  site  and  2,000  which  was  deter¬ 
mined  from  unconsolidated,  undrained  triaxial  shear  test'..  The  undrained  soil 
modulus  was  estimated  using  Eu/su  =  1,000  ,  which  corresponds  to  Eu 
=  2,008,000  psf  or  13,944  psi.  Values  of  Eu  equal  to  3,000,  10,000,  and 
40,000  psi  were  obtained  from  normalized  soil  testing,  pressuremeter  testing, 
and  crosshole  tests.  The  computed  value  of  E^/a^  is  1,050.  Predicted 
f-z  curves  similar  to  the  curves  in  Figure  20  would  be  obtained  using  the 
estimated  soil  properties. 
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Table  1 

Soil-Pile  Adhesion  for  Piles  (After  Tomlinson  1964 


1 


Pile 

Unconfined 

Compressive 

Pile-Soil 

Material 

Strength.  Dsf 

Adhesion,  psf 

Concrete  and 

0  to  1,500 

0  to  700 

timber 

1,500  to  3,000 

700  to  1,000 

3,000  to  6,000 

1,000  to  1,300 

Over  6,000 

Over  1,300 

Steel 

0  to  1,500 

0  to  700 

1,500  to  3,000 

700  to  1,000 

3,000  to  6,000 

1,000  to  1,200 

Over  6,000 

Over  1 , 200 

Table  2 

Comparison  of  Side  Shear  for  UHCC  Teat 


Mea¬ 

sured 

f 

max 

psf 

Computed 

f 

max 

Depth 

ft 

a 

P 

X 

Equa¬ 
tion  20 

Equa¬ 
tion  22 

0 

0 

889 

727 

1,074 

433 

479 

8 

800 

822 

1,100 

1,034 

1 ,008 

1,114 

16 

1  ,400 

862 

1,790 

1,249 

2,110 

2,427 

26 

1 ,250 

667 

1,339 

855 

1,404 

1,535 

30 

1,500 

1,200 

1,380 

2,134 

1 ,641 

1,787 

<13 

Qa  (Computed) 

Q_  (Measured) 

0.83 

1  .  14 

1 .24 

1.22 

1.36 

Table  3 

Comparison  of  Side  Shear  for  SFBM  Test 


Mea¬ 

sured 

Computed 

f 

max 

Depth 

ft 

max 

psf 

q 

8 

X 

Equa¬ 
tion  20 

Equa¬ 
tion  22 

0 

250 

300 

186 

209 

123 

115 

12.5 

250 

325 

344 

333 

246 

229 

17.5 

260 

375 

399 

378 

288 

268 

22.5 

375 

425 

457 

422 

332 

309 

27.5 

475 

475 

520 

466 

379 

353 

32.5 

500 

525 

584 

510 

425 

396 

37.5 

475 

575 

632 

548 

460 

428 

40 

Qs  (Computed) 

Qs  (Measured) 

1.21 

1.18 

1.11 

0.84 

0.79 

Table  4 

Comparison  of  Side  Shear  for  St.  Alban  Teat 


Mea¬ 

sured 

f 

max 


Table  6 

Computed  f-z  Curves 


Method 

Coyle  and 
Reese  (1966) 

Vijay- 

verglya 

(1977b) 

Heydinger  and 
O'Neill  (In 
Preparation 

Kraft,  Ray,  and 
KaRawa  (1981) 

z 

f 

f 

f 

f 

z 

in. 

psf 

psf 

psf 

psf 

in 

APPENDIX  A:  NOTATION 


cm 

*u 

CIU 

D 

E 


f-z  curves 


F 


c 


“o 

ID 

K 


L 
LI 
L L 

■ 


Cross-sectional  area 
Area  of  the  pile  surface 

Cross-sectional  area  of  pile  tip 

Pile-soil  adhesion 

Mean,  undrained  shear  strength 

Average,  undrained  shear  strength 

Consolidated,  isotropic,  undrained  triaxial  testing 
Pile  dianeter 

Young's  modulus  of  pile  resistance 
Undrained  modulus  of  elasticity 

Slope  of  the  initial  tangent  of  curve  f  versus  z/D 
Average  undrained  soil  modulus  for  a  soil  stratum 

Maximum  side  resistance 

Predictive  methods  for  computing  capacity  and  shear  transfer 
versus  pile  movement 

Ratio  of  the  undrained  shear  strength  using  unconsolidated 
undralned  triaxial  tests  to  the  strengths  obtained  by  other 
methods 

Length  correction  factors 

Undrained  shear  modulus 
Initial  shear  modulus 

Secant  shear  modulus 

Inside  diaawter 

Correlation  factor  for  pile  length  to  diameter  ratio  rn 

Hey dinger  and  O'Neill  method 

Ratio  of  the  horizontal  effective  stress  to  the  verr.ca. 
effective  stress  for  the  at-rest  condition 

Pile  length 

Liquidity  index 

Liquid  limit 

Shape  parameter  for  Heydinger  and  O’Neil,  method! 


A  ' 


‘o 


H 

“c 

OD 

OCR 

Pa 

PI 

Ou 

0-2  curve 


Intercept  values  for  m 

(6  sin  ♦' )/(3  -  sin  ♦' ) 

Cohesive  bearing  factor 

Outside  diameter 

Overconsolidated  ratio  of  the  soil 
Atmospheric  pressure 

Plasticity  index 
Total  pile  capacity 

Tip  load-deformation  behavior 

Radial  distance 

Radial  zone  of  influence 


rQ  Radius  of  pile 

Rp  Stress-strain  curve-fitting  constant 

su  Undrained  shear  strength 

(s  \  Undrained  shear  strength  for  soil  at  critical  state 

/cs 

u*  Vertical  pile  movement  resulting  in  maximum  side 
UU  Unconsolidated,  undrained  triaxial  testing 

z  Critical  pile  movement 


a 

e 


( Ao  /u  ) 
rc  e 


°;/o3 


Vertical  displacement  at  the  pile  surface 

Ratio  of  pile-soil  adhesion  to  undrained  shear  strength 

Correlation  factor  for  in  situ  vertical  overburden  pressure  to 
side  resistance  along  a  pile 

Change  in  effective  radial  stress  divided  by  the  excess  pore 
pressure,  u@ 

Parameter  relating  pile  stiffness  to  soil  stiffness 

Mean  effective  overburden  pressure 

Vertical  consolidation  pressure  after  swelling 

Maximum  vertical  consolidation  pressure 

Ratio  of  the  maximum  stress  to  the  minimum  stress 

Shear  stress 

Soil  shear  stress  at  the  pile  surface 


A  2 


Maxi, 


shear  stress  Mobilized  at  pile  failure 


True  angle  of  internal  friction 
Effective  angle  of  internal  friction 

Effective  angle  of  friction  for  soil  at  the  critical  state 
Consolidated,  undrained  triaxial  tests 
Angle  of  friction  between  pile  and  soil 


A3 


